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possibility  of  narrowbeam  detection  of  a  low  frequency  acoustic  signal  wave. 

The  interaction  of  two  plane  waves  is  examined  theoretically  in  a  manner 
similar  to  that  considered  by  Westervelt.  Of  particular  interest  is  the  casi 
in  which  the  pump  wave  is  assumed  to  be  a  plane  wave  originating  at  and  prop¬ 
agating  perpendicular  to  some  planar  boundary.  The  low  frequency  signal  wave 
is  assumed  to  be  present  over  all  space.  The  interaction  of  two  plane  waves 
is  of  interest  for  two  reasons. 

(l)  This  interaction  produces  a  truncated  end- fire  array  (i.e.,  the  parametric 
receiving  array)  with  a  length  equal  to  the  spacing  between  the  planar 
boundary  and  the  observer,  and  (2)  this  interaction  relates  directly  to  the 
problem  of  sound  scattered  by  sound  as  formulated  by  Westervelt. 

Next,  the  solution  for  the  interaction  of  a  high  frequency  spherical  wave 
from  a  point  source  and  a  low  frequency  plane  wave  is  obtained  theoretically 
using  a  two-dimensional  stationary  phase  solution  of  Westervelt' s  scattering 
integral.  The  point  source  solution  is  in  turn  used  to  generate  the  solution 
for  parametric  receiving  arrays  with  various  pump  transducers  including  a 
truncated  line  source  transducer,  a  rectangular  piston  transducer,  and  a 
circular  piston  transducer.  In  each  case,  the  observer  is  assumed  to  be  in 
the  farfield  of  the  pump  transducer.  However,  the  result  includes  the  effects 
of  interaction  in  the  nearfield  of  the  pump  transducer.  The  solution  of  the 
parametric  receiving  array  with  a  point  source  pump  and  a  truncated  line 
receiver  is  also  found.  The  theoretical  solutions  also  include  the  effects 
of  misalignment  of  either  the  pump  or  receiving  transducer.  With  either  the 
pump  transducer  or  the  receiver  misaligned,  the  difference  frequency  beam 
pattern  is  an  asynmetrical  beam  pattern  and  is  the  mirror  image  of  the  sum 
frequency  beam  pattern.  This  property  is  examined  in  detail  in  terms  of  the 
problem  of  two  sound  waves  interacting  at  nonzero  angles. 

A  series  of  experiments  were  conducted  with  a  48  ft  parametric  receiving 
array  with  a  pump  frequency  of  90  kHz.  In  the  experiments,  an  omnidirectional 
transducer  was  used  with  either  a  rectangular  piston  transducer  or  a  square 
piston  transducer.  Each  transducer  could  be  used  as  either  the  pump  trans¬ 
ducer  or  the  receiver  so  that  a  variety  of  parametric  receiving  array  config¬ 
urations  could  be  realized.  A  small  rotator  was  used  to  allow  independent 
rotation  of  the  piston  transducer  for  measurement  of  the  effects  of  misalign¬ 
ment  of  the  piston  transducer.  Theoretical  and  experimental  results  are 
compared  for  the  parametric  receiving  array  with  several  transducer  arrange¬ 
ments.  The  agreement  between  the  theory  and  the  experiment  was  excellent. 
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ABSTRACT 


The  parametric  reception  of  a  low  frequency  plane  wave  by  the 
use  of  nonlinear  interactions  between  acoustic  waves  is  examined  both 
theoretically  and  experimentally.  The  parametric  reception  of  a  low 
frequency  wave  is  accomplished  by  the  use  of  a  high  frequency  acoustic 
pump  wave  which  interacts  with  the  low  frequency  signal  wave  to  produce 
sound  waves  at  the  sum  and  difference  frequencies.  These  sound  waves  are 
received  by  a  second  transducer  placed  on  the  axis  of  the  pump  transducer. 
This  type  of  parametric  array  allows  for  the  possibility  of  narrowbeam 
detection  of  a  low  frequency  acoustic  signal  wave. 

The  interaction  of  two  plane  waves  is  examined  theoretically  in 
a  manner  similar  to  that  considered  by  Westervelt.  Of  particular  interest 
is  the  case  in  which  the  pump  wave  is  assumed  to  be  a  plane  wave  originating 
at  and  propagating  perpendicular  to  some'  planar  boundary.  The  low  frequency 
signal  wave  is  assumed  to  be  present  over  all  space.  The  interaction  of 
two  plane  waves  is  of  interest  for  two  reasons. 

(l)  This  interaction  produces  a  truncated  end-fire  array  (i.e., 
the  parametric  receiving  array)  with  a  length  equal  to  the  spacing  between 
the  planar  boundary  and  the  observer,  and  (2)  this  interaction  relates 
directly  to  the  problem  of  sound  scattered  by  sound  as  formulated  by 
Westervelt. 

Next,  the  solution  for  the  interaction  of  a  high  frequency 
spherical  wave  from  a  point  source  and  a  low  frequency  plane  wave  is 
obtained  theoretically  using  a  two-dimensional  stationary  phase  solution 
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of  Westervelt's  scattering  integral.  The  point  source  solution  is  in  turn 
used  to  generate  the  solution  for  parametric  receiving  arrays  with  various 
pump  transducers  including  a  truncated  line  source  transducer,  a  rectangular 
piston  transducer,  and  a  circular  piston  transducer.  In  each  case,  the 
observer  is  assumed  to  be  in  the  farfield  of  the  pump  transducer.  However, 
the  result  includes  the  effects  of  interaction  in  the  nearfield  of  the  pump 
transducer.  The  solution  of  the  parametric  receiving  array  with  a  point 
source  pump  and  a  truncated  line  receiver  is  also  found.  The  theoretical 
solutions  also  include  the  effects  of  misalignment  of  either  the  pump  or 
receiving  transducer.  With  either  the  pump  transducer  or  the  receiver 
misaligned,  the  difference  frequency  beam  pattern  is  an  asymetrical  beam 
pattern  and  is  the  mirror  image  of  the  sum  frequency  beam  pattern.  This 
property  is  examined  in  detail  in  terms  of  the  problem  of  two  sound  waves 
interacting  at  nonzero  angles. 

A  series  of  experiments  were  conducted  with  a  48  ft  parametric 
receiving  array  with  a  pump  frequency  of  90  kHz-  In  the  experiments, 
an  omnidirectional  transducer  was  used  with  either  a  rectangular  piston 
transducer  or  a  square  piston  transducer.  Each  transducer  could  be  used 
as  either  the  pump  transducer  or  the  receiver  so  that  a  variety  of  parametric 
receiving  array  configurations  could  be  realized.  A  small  rotator  was  used 
to  allow  Independent  rotation  of  the  piston  transducer  for  measurement  of 
the  effects  of  misalignment  of  the  piston  transducer.  Theoretical  and 
experimental  results  are  compared  for  the  parametric  receiving  array  with 
several  transducer  arrangements.  The  agreement  between  the  theory  and  the 
experiment  was  excellent. 
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I.  INTRODUCTION 


It  is  veil  known  that  a  finite  amplitude  sound  vave  will  distort 

as  it  propagates.  "  This  nonlinear  propagation  or  self- interaction 

occurs  because  the  phase  velocity  is  a  function  of  the  instantaneous 

amplitude  of  the  particle  velocity.  It  would  also  be  expected  that  two 

sound  waves  of  frequencies  and  fg  present  simultaneously  in  a  medium 

would  generate  sum  and  difference  frequencies,  f  ±  fg.  A  spurious 

difference  frequency  tone  was  observed  by  a  German  organist,  Sorge,  in  1745, 

and  later  in  1754  by  an  Italian  violinist,  Tart  ini.  ^  These  Bpurious 

frequency  components  have  been  known  to  musicians  ever  since  and  are 

referred  to  as  "combination  tones."  In  1875*  Helmholtz1®  predicted  that  a 

sum  frequency  tone  should  also  exist.  With  the  ubb  of  a  Helmholtz  resonator, 

Helmholtz  showed  that  the  difference  frequency  tone  was  generated  by  such 

instruments  as  the  harmonium  when  two  high  frequency  tones  were  played 

simultaneously.  Helmholtz  suggested  that  this  tone  was  generated  in  the 

air  within  the  instrument.  On  the  other  hand,  the  presence  of  two  notes 

from  a  pipe  organ  had  little  effect  on  a  resonator.  In  this  case,  Helmholtz 

concluded  that  the  difference  frequency  was  either  subjective  or  generated 

within  the  ear.  The  existence  of  the  sum  frequency  was  a  matter  of  debate 

until  1895  when  Rucker  and  Edser11  observed  its  effect  on  a  tuning  fork 

12 

of  the  same  frequency.  Rayleigh  discussed  the  subject  in  his  book 

15  2 

The  Theory  of  Sound  in  1910.  Lamb  used  Airy's  method  of  successive 
approximations  to  demonstrate  the  existence  of  sum  and  difference 
frequency  components.  The  existence  of  sum  and  difference  frequency 
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lk 

components  in  the  throat  of  a  horn  was  considered  by  Rocard  in  1933  and 

4 

by  Thuras,  Jenkins,  and  O'Neil  in  1935*  Thuras,  Jenkins,  and  O'Neil 
measured  the  sum  and  difference  frequency  components  experimentally  using 
a  long  tube. 

The  interaction  of  two  sound  waves  in  a  bounded  region  was 

considered  again  in  1956  by  Ingard  and  Pridmore-Brown.^  In  1957, 

17 

Westervelt  *  generated  a  solution  for  the  interaction  of  two  infinite 

plane  waves  with  no  boundaries  present. 

Despite  this  activity,  the  interaction  of  sound  waves  did  not 

l6 

receive  widespread  attention  until  the  early  1960's  when  Westervelt 
proposed  the  nonlinear  or  parametric  acoustic  array  for  the  transmission 
or  detection  of  a  low  frequency  acoustic  wave.  He  showed  that  a  highly 
directional,  low  frequency  sound  beam  could  be  generated  with  a  Bmall 
transducer  by  the  use  of  the  interaction  of  two  high  frequency  carrier 

sound  beams.  This  parametric  transmitting  array  has  been  demonstrated 

19  20  21 

experimentally  by  Beilin  and  Beyer  ,  Hobaek,  Zverev  and  Kalachev, 

22  25  24  25 

Muir  and  Blue,  Smith,  Merklinger,  Truchard  and  Wlllette,  among 


others. 


26-28 


The  parametric  reception  of  a  low  frequency  wave,  also  proposed 
18 

by  Westervelt,  is  accomplished  by  the  use  of  a  high  frequency  acoustic 

pump  wave  which  interacts  with  the  low  frequency  signal  wave  to  produce 

sound  waves  at  the  sum  and  difference  frequencies.  These  sound  waves  are 

received  by  a  second  transducer  placed  on  the  axis  of  the  pump  transducer. 

29  30 

The  parametric  receiving  array  was  first  studied  by  Berktay  and  later 
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by  Tjotta.  An  experiment  was  conducted  by  Berktay  and  Al-Temimi7 

that  demonstrated  the  existence  of  the  parametric  receiving  array  and, 


! 


3 


furthermore  ,  shoved  that  the  array  had  a  directivity  function  which  is 

very  similar  to  that  of  a  conventional  end-fire;  array  with  a  length 

equal  to  the  distance  between  the  two  transducers.  Similar  experiments 

34  35 

have  been  conducted  by  Barnard  et  al.,  Konrad  et  al.,  Muir  and 
Berktay,*^  Zverev  and  Kalachev, ^  and  Date  and  Tozuka.^ 

Most  of  the  experimenters  developed  a  theoretical  model 
designed  to  correspond  to  the  conditions  of  the  particular  experiment. 

The  experiments  conducted  by  Berktay  and  Al-Temimi  were  primarily  done  in 
the  nearfield  of  the  transducer.  Consequently,  the  theory  was  appropri¬ 
ately  designed  for  consideration  of  interaction  in  the  nearfleld  of  the 
pump  transducer.  On  the  other  hand,  the  numerical  solution  by  Barnard 
et  al.  assumed  that  the  interaction  occurred  only  in  the  farfield  of  the 
pump  transducer.  The  theory  used  by  Muir  and  Berktay  for  the  analysis 

of  a  single  parametric  array  is  described  fully  in  a  paper  by  Berktay  and 
30 

Shooter.  A  closed  form  solution  vas  found  for  the  sum  and  difference 

frequency  sound  field  for  the  parametric  receiving  array  in  the  farfield 

of  the  narrcwbeam  pump  transducer.  The  theoretical  approach  of  Zverev  and 

Kalachev  and  of  Date  and  Tozuka  does  not  use  Westervelt's  perturbation 

procedure  for  finding  the  value  of  the  sum  and  difference  frequency  sound 

pressure.  Instead,  the  modulation  of  the  high  frequency  pump  wave  by  the 

low  frequency  signal  wave  is  considered  in  terms  of  the  phase  modulation 

of  the  pump  wave  caused  by  the  time  varying  change  in  the  sound  speed  due 

4o 

to  the  presence  of  the  signal  wave.  Recently,  Rogers  et  al.  obtained 
a  solution  for  the  parametric  receiving  array  in  the  nearfield  of  a  small 
narrowbeam  pump  transducer. 


4 


In  each  of  the  models  described  above,  the  effects  of  shadowing 

of  the  signal  wave  by  the  pump  transducer  were  neglected.  These  effects 

4l 

were  considered  by  Al-Temimi  and  were  found  to  reduce  the  directivity 

of  the  parametric  receiving  array  when  the  pump  transducer  was  large 

compared  to  the  signal  frequency  wavelength. 

The  possibility  for  utilizing  an  arrangement  similar  to  the 

parametric  receiving  array  for  the  construction  of  a  traveling  wave 
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parametric  amplifier  was  proposed  by  Tucker,  Stepanov,  and  Berktay. 

However,  each  researcher  concluded  that  amplification  of  the  signal  wave 

is  not  possible  unless  the  sum  frequency  wave  is  suppressed  or  the  medium 

is  dispersive.  Such  amplification  in  a  dispersive  waveguide  is  predicted 

44 

by  Ostrovskii  and  Papilova. 

45 

Berktay  and  Al-Temimi  have  recently  related  the  parametric 

receiving  array  to  the  problem  of  sound  scattered  by  sound.  This 

relationship  is  an  interesting  one  and  deserves  attention  because  it  has 

continually  been  debated  since  Ingard  and  Pridmore-Brcwn  first  studied 
15 

the  problem.  Ingard  and  Pridmore -Brown  predicted  that  a  scattered 
sound  field  would  propagate  outside  the  interaction  region  common  to  the 
two  primary  waves.  Their  experimental  measurement  seemed  to  verify  the 
predictions.  However,  later  attempts  at  verifying  these  predictions  have 
proved  unsuccessful.  ^,47,48  jn  ^7,  Westervelt^'1^  considered  the 
nonlinear  interaction  of  two  infinite  plane  waves  intersecting  at  nonzero 
angles  and  found  that  the  scattered  pressure  was  related  only  to  the 
value  of  the  primary  sound  wave  amplitude  at  the  observer  point.  This 
solution  did  not  include  the  effects  of  boundaries  or  absorption.  Conse¬ 
quently,  the  solution  had  a  singularity  when  the  primary  sound  waves  were 


5 


propagating  in  the  same  direction.  This  result  was  consistent  with 

12  1^ 

previous  predictions  by  Lord  Rayleigh  and  Lamb  ^  which  had  been  verified 
experimentally  by  Thuras,  Jenkins,  and  O'Neil,  provided  we  assume  that 
the  interaction  occurs  over  an  infinite  distance.  Furthermore,  Westervelt 
predicted  that  no  scattered  sound  exists  outside  the  region  of  interaction. 

Various  magnitudes  of  scattered  sound  pressure  outside  the  region  of 

4q  .  50 

interaction  have  been  predicted  by  Dean,  Lauvstad  and  TJotta,  and 

47 

Al-Temimi.  ' 

In  1972, Westervelt  recast  the  1957  result  into  a  form  which 
permitted  the  singularity  to  be  removed.  He  furthermore  related  this 
new  solution  to  the  truncated  end- fire  array,  i.e.,  the  parametric 
receiving  array  for  the  case  when  the  frequency  of  one  wave  was  considerably 
higher  than  the  other. 

The  present  study  includes  theoretical  analysis  and  experimental 
investigation  of  the  parametric  receiving  array.  The  investigation  of  the 
parametric  receiving  array  is  extended  to  include  omnidirectional  pump 
transducers  and  pump  transducers  which  have  one  dimension  up  to  one  wave¬ 
length  or  more  at  the  signal  frequency  and  two  dimensions  which  are  small 
compared  to  the  signal  frequency  wavelength  in  order  to  minimize  the 
effects  of  acoustic  shadowing.  A  closed  form  solution  of  Westervelt's 
scattering  integral  for  an  omnidirectional  pump  transducer  will  be  found 
and,  in  turn,  used  to  obtain  the  solution  for  a  line  pump  transducer.  The 
theory  is  based  on  Westervelt 1  s  quasilinear  solution  in  which  only  second- 
order  interactions  are  considered.  The  effects  of  absorption  are  ignored 
in  the  derivation  of  the  basic  second-order  wave  equation.  However,  the 
attenuation  of  the  first-  and  second-order  sound  waves  due  to  absorption 
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Is  Included  In  the  description  of  these  sound  waves  In  a  manner  similar 
to  that  used  by  Westervelt.  Westervelt's  solution  for  the  problem 
of  sound  scattered  by  sound  is  analyzed  and  related  to  the  parametric 
receiving  array  for  two  plane  waves.  An  end- fire  array  function  is 
generated  from  this  solution.  The  parametric  receiving  array  offers  an 
excellent  opportunity  to  study  the  scattering  of  sound  by  sound  since 
interaction  occurs  at  nonzero  angles  between  the  direction  of  propagation 
of  the  two  waves.  In  a  manner  similar  to  that  used  by  Westervelt  in  1972, 
the  singularity  is  removed  from  the  1957  result  by  the  addition  of  a 
solution  of  the  homogeneous  wave  equation,  thereby  placing  a  boundary 
condition  on  the  second-order  pressure.  The  second-order  sound  field  is 
shown  to  be  nonzero  even  when  the  two  sound  waveB  are  not  propagating  in 
the  same  direction.  In  Chapter  III,  the  significant  features  of  this 
solution  are  interpreted  in  terms  of  the  properties  of  the  parametric 
receiving  array. 

Also  in  Chapter  III,  the  solution  for  the  interaction  of  a 
spherical  wave  from  a  point  source  and  a  plane  wave  is  found  using 
Westervelt's  scattering  integral.  The  second-order  sound  pressure  is 
found  by  integrating  the  freefield  Green's  function  over  the  volume 
distribution  of  sources  to  form  the  scattering  integral.  This  scattering 
integral  is  in  turn  used  to  find  the  second-order  sum  and  difference 
frequency  sound  pressure.  The  integral  is  evaluated  using  a  two-dimensional 
stationary  phase  integral  solution  technique  which  has  been  used  to  solve 
certain  optics  problems.  This  solution  is  found  to  have  the  same  end-fire 
array  properties  as  the  solution  for  the  interaction  of  two  plane  waves. 
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In  Chapter  IV,  the  point  source  solution  is  used  as  a  starting 
point  to  solve  for  the  second-order  sound  field  generated  by  a  truncated 
line  source  and  a  plane  wave;  the  line  source  was  chosen  because  the 
experimental  and  theoretical  models  can  be  made  to  coincide  most  easily. 

The  effects  of  acoustic  shadowing  of  the  low  frequency  plane  wave  are 
minimized  even  when  the  pump  has  one  dimension  on  the  order  of  a  wave¬ 
length  at  the  signal  frequency.  The  large  dimension  for  the  pump  allows 
the  study  of  the  effects  of  pump  size  on  array  beamwidth  and  sensitivity. 
Using  the  large  pump,  the  scattering  of  sound  by  sound  can  also  be  investi¬ 
gated. 

Lastly,  Chapter  IV  considers  the  interaction  of  plane  sound 
waves  and  waves  produced  by  a  pump  transducer  that  is  small  in  all  three 
dimensions  compared  to  a  wavelength  at  the  signal  frequency.  Again,  the 
point  source  solution  is  used. 

In  Chapter  V,  the  experimental  arrangement  is  described.  The 
experiments  were  conducted  aboard  the  STEP  Barge  at  ARL's  Lake  Travis 
Test  Station.  The  geometry  for  the  theoretical  and  experimental  studies 
was  such  that  the  results  could  be  readily  compared.  Several  types  of 
transducers  were  used  to  reproduce  the  geometries  prescribed  by  the 
theoretical  models.  The  arrangement  of  electronic  equipment  is  essentially 
the  same  as  that  used  by  the  author  and  described  in  a  paper  by 

54 

Barnard  et  al. 

Chapter  VI  includes  a  comparison  of  experimental  and  theoretical 
results.  Beam  patterns  with  the  pump  aligned  and  misaligned  are  compared 
in  the  same  coordinate  systems.  Sound  pressure  level  predictions  are 
compared  with  the  experimental  results.  The  agreement  between  the  theory 
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and  the  experiment  was  excellent.  The  results  are  discussed  in  light  of 
the  problem  of  sound  scattered  by  sound. 

In  the  final  chapter,  a  summary  of  the  theoretical  and 
experimental  results  of  the  present  study  is  given. 


II.  THE  EQUATIONS  OF  HYDRODYNAMICS  FOR  THREE-DIMENSIONAL 
WAVE  MDTION  IN  A  LOSSLESS  FLUID 

In  elementary  treatments  of  acoustic  wave  propagation  it  is 
usually  assumed  that  waves  of  infinitesimal  amplitude  are  being  modeled 
so  that  the  intrinsically  nonlinear  acoustic  equations  can  be  linearized. 
Using  this  assumption,  acoustic  wave  propagation  without  distortion  is 
predicted.  Likewise,  if  acoustic  waves  with  two  or  more  frequencies  are 
present,  no  intermodulation  is  predicted  for  the  propagation  of  these 
waves.  If  the  nonlinearity  of  the  acoustic  equations  is  taken  into 
account,  an  originally  sinusoidal  waveform  will  indeed  distort  as  it 
propagates  or,  if  two  waves  are  present,  intermodulation  products  will  be 
predicted.  These  intermodulation  products  will  be  present  whether  the 
amplitude  of  the  two  waves  is  large  or  small.  The  ability  to  detect  these 
lntermcdulation  products  will  depend  on  whether  the  amplitude  is  large 
enough  to  be  detected  in  the  presence  of  the  background  noise.  In  order 
to  account  for  this  intermodulation,  the  elementary  wave  equation  must  be 
modified  to  take  into  account  higher  order  terms  present  in  the  propagation 
of  a  wave. 
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A  perturbation  method  used  by  Eckart,  and  later  by  Lighthill  '^ 

XT 

and  Westervelt,  *  can  be  used  to  obtain  a  wave  equation  which  includes 

second-order  components  to  provide  a  better  model  for  the  propagation  of 

a  finite  amplitude  sound  wave.  In  this  theory,  sometimes  referred  to  as 

small- signal  nonlinear  theory  or  quasi  linear  theory,  only  the  first 

modulation  product  terms  are  considered  with  a  second-order  wave  equation 

which  has  source  terms  proportional  to  the  time  derivative  of  the  square 
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of  the  first-order  sound  pressure.  If  two  sound  fields  with  two  different 
frequencies  are  present,  a  source  term  is  predicted  that  is  proportional 
to  the  time  derivative  of  the  cross  product  of  the  two  first-order  sound 
pressures.  In  this  manner,  either  of  the  two  sound  fields  can  be  a 
superposition  of  several  waves  each  having  the  same  frequency.  It  must 
be  noted  that  this  model  is  valid  only  if  the  two  waves  have  amplitudes 
small  enough  that  the  higher  order  terms  can  be  ignored.  In  other  words, 
if  excess  attenuation  or  amplitude  loss  is  present  due  to  finite  amplitude 
effects  this  model  will  no  longer  be  valid.  However,  a  large  number  of 
problems  concerning  the  parametric  receiving  array  can  be  solved  by  using 
the  second-order  or  quasilinear  solution. 

The  properties  of  a  fluid  can  be  described  in  either  Eulerian 

or  Lagrangian  coordinates.  The  relationship  of  these  two  coordinate 
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systems  with  real  sources  present  has  been  studied  by  Kline.  In  the 
present  study,  the  Eulerian  coordinate  system  1b  used  throughout. 

Whenever  necessary,  the  boundary  conditions  will  be  simplified  so  that 
the  boundary  conditions  can  be  expressed  in  Eulerian  coordinates  without 
difficulty. 

We  shall  derive  a  second-order  wave  equation  in  a  manner  similar 
57 

to  that  used  by  Westervelt.  The  equations  for  the  conservation  of  mass 
and  momentum  will  be  given  with  real  sources  and  viscosity  included. 
However,  absorption  and  real  sources  will  be  ignored  in  the  derivation 
of  the  wave  equation.  The  effects  of  absorption  will  be  included  in  an 
"ad  hoc"  manner  in  the  solution  of  the  second-order  wave  equation  in 
Chapter  III.  If  the  effects  of  real  sources  are  to  be  included,  the  first- 
order  substitution  into  the  second-order  wave  equation  must  include  these 
sources.  The  equation  of  state  will  be  that  of  an  isentropic  fluid. 
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In  Eulerian  coordinates,  the  equations  of  the  conservation  of 

58 

mass  and  momentum  with  sources  present  can  he  expressed  as  follows: 

Conservation  of  Mass 

dp, 


T 

“5t  + 


^•PTU  = 


(2.1) 


Conservation  of  Momentum 


dpTu 

~3T 


+  V*T 


=  t 


(2.2) 


where  is  the  total  density,  u  is  the  particle  velocity,  pQ  is  the  ambient 
or  undisturbed  density,  Q  is  the  rate  of  introduction  of  new  fluid  mass  into 
a  unit  volume,  T  is,  a  tensor  such  that  T  is  the  sum  of  Tg,  the  stress 
tensor,  and  V  the  momentum  flux  tensor,  and  F  1b  the  external  force 
per  unit  volume. 

The  momentum  flux  tensor  can  be  expressed  as 


/pTuxUx 

pTUxUy 

pTUxUi 

pu  u 

y  x 

pu  u 
y  y 

pu  u 
y  z 

\puzux 

pu  u 
z  y 

pu  u 
z  z 

where  u  ,  u  ,  and  u  are  the  three  components  of  u. 
x  y  z 

The  stress  tensor  can  be  expressed  as 


(2.5) 


(2.4) 


where  p  is  the  total  pressure  and 


I 


»u  -  -  f)  « -  a.  ^ 

Bu  -  +  S*) 
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(2.5) 


i  *  J 


where  p  is  the  coefficient  of  shear  viscosity  and  rj  is  the  coefficient  of 
bulk  viscosity. 

For  inviscid  fluids,  the  tensor  Tg  is  a  simple  one  with  p  on  the 
main  diagonal  and  zeroes  elsewhere.  In  this  case,  ^»Tg  is  simply  Vp.  To 
obtain  the  expression  for  we  consider  the  x  component  such  that 


a(pTUx}  a(pTUv)  a(pTuz)  du- 

(  V.  Tm)  =»  U  '  1  ■  +  u  -  \  ■  +  U  — »■  +  P  u 

v  M'x  x  ox  x  py  x  dz  T  : 


x 

X  dx 


du  du 

+  PTUy  cbc^  +  PTUz  =  \  (pT^  +  PT^’^  ux  * 


(2.6) 


Then  we  have 


^•Tm  =  uv*  (ptu)  +  pT(u*V)u  . 


(2.7) 


The  momentum  equation  for  an  inviscid  fluid  now  can  be  expressed  as 


dpTu 

“5 r 


r-  +  $P  +  u  &  (p-u)  +  pju.^u  =  t 


(2.8) 


The  acoustic  pressure  can  be  related  to  the  density  through  the 
equation  of  state: 

p  =  F(pT,S)  ,  (2.9) 

where  S  is  the  entropy. 

For  an  isentropic  fluid  we  have 

P  =  f(pt) 


(2.10) 
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We  shall  use  a  perturbation  procedure  similar  to  that  used  by 
Eckart,*^2  Lighthill, ^  and  Westervelt.^*^  The  field  variables  u)  p^,  and 

p  can  be  expanded  in  series  such  that: 

—*  .  2— > 

U  =  \u^  +  X  Ug  +  .  .  .  , 

PT  =  P0  +  Xp^^  +  X2P2  +  •  •  •  t  and  (2  •  n) 

2 

p  =  Po  t  \p1  +  X  p2, 

where  Ms  a  nondimensional  mathematical  parameter  introduced  for  con¬ 
venience  and  p  is  the  ambient  pressure.  We  have  assumed  that  u  =0.  An 
o  o 

approximate  equation  of  state  for  an  isentropic  fluid  can  be  obtained  by 
expanding  p  in  a  Taylor  series  and  keeping  only  the  first  two  terms  such 
that 


=  p° +  fel  0(pT'p°)+  ©  (Pt'p°>! 

S,°  T  B,0 


where  S  is  the  entropy. 


The  infinitesimal  sound  speed  cq  is  defined  by 


(2.32) 


(2.13) 


The  pressure  can  also  be  expressed  as 

p  =  p0  +  cQ2  (PT-PQ)  +  \  (PT-P0)2  , 


(2.1*0 


where 


B/A  =  P0/co2(^  p/3pT2) 


S,0 


(2.15) 


By  substituting  Equation  (2.1l)  into  Equation  (2.l4)  and  assembling  like 


powers  of  \,  we  get 


pl  =  co  pl 


(2.16) 


2  J 
c  p.  + 
o  y2 


=  co%  + 


(b_\  fo_  2 

Wp0  Pl 

(pi} 

V2A/  „  2 

\  /  o  c 


(2.17) 


If  we  assume  that  the  source  terms  in  Equations  (2.1)  and  (2.8)  are  first- 
order  functions,  we  can  obtain  the  first-order  equations  by  substituting 
Equations  (2.1l)  into  (2.1)  and  (2.8),  keeping  terms  with  first  powers 
of  \.  We  have  the  following  equations: 

Conservation  of  Mass 


Spl 

"St;  +  Po^*Ul  =  ^ 


(2.18) 


Conservation  of  Momentum 

do, 

po“5t  +  \  =  ^ 

Wave  Equation  for  p^,  p^,  aid  u 


(2.19) 


--  -  c  =  p  &  - 

.2  o  H1  Ko  dt 


or  ^p  -  - - 1  =D2p1  =  -P  ^  , 

1  c  2  dt2  1  0  ^ 


where  □ 


.j.jL-i. 


c  2  dt2 
o 


is  the  wave  operator. 


(2.20) 
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and 


„2->  1  df  .  & 

Qi^  =  -  - g  ^  +  T3 


P  c 
po  o 


Likewise,  the  second-order  equations  can  be  obtained  by  equating 
2 

terms  containing  \  to  obtain  the  following  second  order  equation. 

Conservation  of  Mass 


^2 

It  +  PoV,u2  +  V*plul  =  0 


(2.21) 


Conservation  of  Momentum 


po  ""5t  +  +  PA  (^*^1)  +  p0  (V^i  +  %  =  0  (2,22) 


Wave  Equation 

=O^P2  -  ^  +  P0(V^l] 

+  v2(P2-eo2p2) 


(2.23) 


If  we  assume  i?  is  irrotational,  we  can  find  a  scalar  0  such  that 
we  have  n=^>.  We  shall  denote  single  time  derivatives  with  a  dot  above 
the  parameter  and  double  time  derivatives  with  two  dots.  If  we  furthermore 
assume  that  no  real  sources  are  present,  the  second  order  wave  equation 
can  be  expressed  as 
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ifpa  -  -p0[(*^1)2  +  (%)-j  (%>  .  ±  v5  (%•%)] 

^(Pl>2 
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=  -p, 


!-§- T  (p/  +  +  lei2  (%•%) 

Lp0  c0 


_1_£_ 

2c  2  dt2 
o 


1  IB' 

iJT'W' 

0  o 


(^%)j 

^(px)2 


(2.25) 


dt‘ 


We  can  use  the  following  identities  to  simplify  the  wave  equation: 


,  a2  L)\ 

I  (%•%)='  -^r ♦  <%)•  (^r-) 


•  ^Pl  •  ^P,  - 


(2.26) 


(%  (V2^x)  =  -%  ^  .  — £±  = 


1  T*A 


^  (%>  _ 


dtc 


P  C 
Mo  o 


•  ^P-L 


(2.27) 


(p0x)  =  2(Vpx)  •  V>01  - 


(2.28) 


□2pl2  =  — ~ri  (Pi^2  +  2^!*^! 


(2.29) 


and 
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Then  the  wave  equation  for  is: 


A  .  i  h  \  tlhl 

u  P«  -  5:  ^2Aj  2 


or 


P  c 
Ko  o 


dt‘ 


-P, 


r  1  , 1  (1  =,  * 

“a~f  <I1)  +  TT l~a  - 5--) 

Lpo  °o  co  Vo  0  / 

-4 


P  c 

°°  -  v  g(  2n 

1  /b  \  r(pi  > 
T  \2A ) . Z2 


p  c 
o  o 


dt£ 


2  ,2  .  2  ^ 1 > 


-pol7T-&(V  +77T  ^2 

Ko  o 


p2  =  - 


[2 

(Z4._Z4.!°%%Y 

V  Co  2poCo  / . 

1  A  B  \  ^pl2 

7M1  +  ^1F 


O  O 


^  r  1  a  •  1  21  ->-> 

U  I  — {O^p^ - 2  P1  "  2  Po  ul’ul 

I  c  2p  c 

L  o  Ko  0  J 


The  second-order  wave  equation  can  alternatively  be  expressed  in  the 
following  form: 


p2  =  'po 


c 

L  0 


1 

2 


+ 


dtc 


2  "  po  2 


vv  ** 

o  0 


"u 

1  1 


.  2p  c 
Loo 


(2.30) 


(2.31) 


(2.32) 
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The  first  two  terms  on  the  right-hand  side  of  Equations  (2 . 3l)  and  (2.32) 
are  the  source  terms  for  the  second-order  pressure.  The  first  term  in  each 
expression  is  present  whether  or  not  the  fluid  is  nonlinear.  The  second 
term  is  due  to  the  nonlinearity  of  the  fluid.  The  term  • )  on  the  right- 
hand  side  of  Equations  (2.31)  and  (2.32)  is  also  present  whether  or  not  the 
fluid  is  nonlinear.  The  term  □2(  •  )  is  ignored  by  Westervelt^  in  the 
derivation  of  the  second-order  wave  equation  since  the  term  yields  a  local 
pressure.  Therefore,  it  is  not  cumulative  and  dots  not  propagate  outside 
the  interaction  region.  For  the  parametric  transmitter,  the  low  frequency 
sound  wave  generated  by  the  interaction  of  two  high  frequency  waves  will 
propagate  outside  the  region  of  interaction  and  the  term  can  usually  be 
ignored.  For  the  parametric  receiving  array,  the  observer  is  in  the 
region  of  interaction;  therefore,  the  term  can  only  be  ignored  if  it  is 
small  compared  to  the  first  two  terms.  The  physical  significance  of 
Equation  (2.32)  can  be  more  readily  understood  if  we  obtain  an  expression 
for  the  second-order  pressure  in  terms  of  the  second-order  scalar  poten¬ 
tial  We  can  obtain  such  an  expression  by  integrating  the  second- 
order  momentum  equation  over  space  and  using  the  appropriate  first-order 
relations  to  obtain: 


^2 

p2  *  •  p0 


where 


1  /  -  -  pl  \ 

"  2  lpoVul - 2  ) 

\  PoCo  / 


(2.33) 


is  a  local  pressure  due  to  the  presence  of  the 


first-order  sound  field.  The  second  wave  equation  (Equation  2.32)  includes 


this  pressure  in  a  simple  form. 
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In  Chapter  III  we  shall  solve  the  second-order  wave  equation 
(Equation  2.31  or  Equation  2.32)  for  several  configurations  of  the  first- 
order  sound  field.  The  first  of  these  configurations  will  be  two  infinite 
plane  waves  and  then  two  plane  waves  present  over  a  half  space.  We  shall 
show  that  Equations  (2.31)  and  (2.32)  are  equivalent  for  a  first-order 
field  consisting  of  two  plane  waves.  We  shall  also  use  a  Green's  function 
solution  of  Equation  (2.3l)  to  solve  for  the  second-order  pressure  gener¬ 
ated  by  the  interaction  of  a  spherically  spreading  wave  and  a  plane  wave. 
For  this  solution  we  shall  only  consider  the  source  term, 


<55 

Westervelt . 


in  a  manner  similar  to  that  suggested  by 


III.  SOLUTION  OF  THE  SECOND-ORDER  WAVE  EQUATION 


In  the  last  chapter,  we  derived  the  second-order  wave  equation 
with  a  source  function  which  is  quadratic  in  the  first-order  field 
variables.  For  the  parametric  receiving  array,  we  are  interested  in  the 
sum  and  difference  frequency  sound  pressures  which  are  generated  by  the 
cross-product  terms  in  the  source  function. 

In  Section  A,  we  present  a  simple  means  to  generate  the  source 
function  for  the  sum  and  difference  sound  pressures  by  the  use  of  complex 
numbers. 

In  Section  B,  we  consider  the  interaction  of  two  plane  waves 
with  cm  angle  6  between  their  directions  of  propagation.  The  interaction 
of  two  infinite  plane  waves  is  of  interest  because  the  problem  relates 
directly  to  the  problem  of  the  scattering  of  sound  by  sound  aB  first 
formulated  by  Westervelt.  For  the  lossless  case,  the  source  function  is 
equal  to  the  D'Alembertian  of  a  function  of  the  first-order  sound  fields. 
Hence,  the  value  of  the  second- order  sound  pressure  is  dependent  only  on 
the  local  sound  pressure  unless  the  two  first-order  sound  waves  are 
propagating  in  the  same  direction.  In  the  case  of  collinear  waves, 
the  predicted  value  for  the  sound  pressure  for  the  sum  and  difference 
frequency  components  is  infinite  as  would  be  expected  if  no  boundaries, 
absorption,  or  finite  amplitude  effects  are  present.  This  result  led 

Westervelt  to  predict  that  no  scattering  of  sound  by  sound  would  exist  if 

\ 

the  two  sound  waves  were  not  propagating  in  the  same  direction.  In  the 
last  part  of  Section  B,  we  shall  consider  the  interaction  of  one  sound 
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field  present  over  all  space  and  a  second  present  over  one-half  space 
starting  at  some  planai  boundary  and  propagating  away  from  and  perpen¬ 
dicular  to  this  boundary.  To  solve  this  problem,  we  add  a  homogeneous 
solution  of  the  wave  equation  to  the  solution  obtained  in  the  first  part 
of  Section  B.  In  this  derivation,  we  shall  neglect  the  effects  of  real 
sources  that  would  be  necessary  to  generate  the  first-order  fields.  In 
fact,  we  shall  place  the  second-order  pressure  equal  to  the  local  pressure 
expressed  in  Equation  (2.33)  at  the  planar  boundary.  If  we  observe  the 
second-order  pressure  at  some  distance  from  the  boundary,  we  find  that 
this  pressure  is  cumulative  and  propagates  in  a  direction  different  from 
the  direction  of  either  of  the  primary  waves.  The  dependence  of  the 
pressure  amplitude  on  the  angle  between  the  propagation  vectors  of  the 
two  primary  waves  is  very  similar  to  the  directivity  function  of  the 
parametric  receiving  array  as  found  by  Berktay. 

For  two  plane  waves  we  have  no  difficulty  in  solving  the  wave 
equation  directly.  However,  in  general,  we  will  find  it  necessary  to 
formulate  the  solution  in  terms  of  an  integral  of  a  Green's  function.  A 
description  of  the  Green's  function  solution  is  given  in  Section  C. 

In  order  to  solve  the  problem  of  a  parametric  receiving  array 
with  an  arbitrarily  shaped  planar  piston  pump,  we  shall  find  a  solution 
which  can  be  used  as  a  starting  point  for  solving  more  complex  problems. 
The  most  obvious  solution  is  the  solution  of  a  hi^h  frequency  point  source 
pump  and  a  plane  wave  low  frequency  signal.  If  we  have  the  solution  for 
an  omnidirectional  pump  transducer,  the  solution  for  virtually  any  other 
type  of  two-  or  three-dimensional  pump  transducer  can  be  obtained  by  sum¬ 
ming  this  solution  over  the  active  face  or  volume  of  the  transducer. 
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It  is  in  this  light  that  we  shall  generate  the  solution  for  the 
Interaction  of  a  spherical  wave  from  a  point  source  and  a  plane  wave.  We 
shall  use  the  Green's  function  integral  solution  for  the  second-order  wave 
equation  with  a  first-order  sound  field  consisting  of  a  spherical  source 
and  a  plane  wave.  This  solution  yields  a  rather  unwieldy  integral  which 
does  not  succumb  to  simple  integration  techniques.  This  integral  can  be 
solved  using  a  two-dimensional  stationary  phase  technique  similar  to  one 
used  in  optics.  We  shall  ignore  the  effects  of  real  sources  and  the  local 
pressure.  In  reality,  the  second-order  pressure  contribution  at  the 
source  and  receiver  may  not  be  negligible  for  a  particular  application  of 
the  parametric  receiving  array  since  we  may  desire  to  use  the  properties 
of  the  directivity  function,  particularly  in  the  nearfield  of  the  parametric 
array,  to  eliminate  an  undesired  noise  source.  If  an  additional  pressure 
term  is  present  due  to  the  pump  source  or  the  local  pressure,  the  perform¬ 
ance  of  the  parametric  receiving  array  may  be  degraded  for  an  application 
requiring  a  very  high  rejection  of  a  target  at  an  undesired  bearing.  In 
order  to  analyze  the  effects  of  the  real  sources  and  the  local  pressure 
terms,  a  thorough  analysis  of  the  particular  type  of  transducer  is  needed. 
Since  the  transducers  are  designed  to  operate  at  the  carrier  frequency, 
the  transducer's  properties  at  the  low  frequency  may  be  such  that  the 
transducer  actually  moves  with  the  low  frequency  fluid  particle  displace¬ 
ment.  This  condition  could  significantly  change  the  way  the  parametric 
array  would  perform  as  compared  to  its  performance  under  the  assumption 
that  the  transducer's  position  remains  fixed.  In  this  study,  this  particular 
problem  will  not  be  considered;  however,  it  does  appear  to  be  worthy  of 


further  investigation. 
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A.  Use  of  Complex  Numbers  ( 

In  Equation  (2.31),  we  can  see  that  the  sum  and  difference 
frequency  source  terms  are  generated  by  cross-products  of  each  of  the  two 
first-order  sound  fields.  It  Is  convenient  to  use  complex  numbers  to 
represent  the  oscillatory  part  of  the  various  functions.  The  real 
physical  variable  is  equated  to  the  real  part  of  the  complex  variable. 

In  the  case  where  quadratic  terms  are  encountered,  special  care  must  be 
taken  in  their  calculation.  For  two  complex  numbers  M  and  N 

Re[M].Re[N]  =  |  Re  M[N+N*] 

where  Re[M]  represents  the  real  part  of  the  complex  variable  and  the 
asterisk  denotes  the  complex  conjugate.  Hie  transformation 

MN  ^  i  MN  +  i  MN* 

2  2 

will  be  used  to  generate  the  quadratic  terms  when  complex  variables  are 
used. 

For  example,  consider  the  case  of  two  plane  waves  with  wave 
vectors  intersecting  at  angle  0.  The  first-order  pressure  is 


pl  "  Plle 


j(kx»  r-cUjt)  j(l^.r-a^t) 

+  P! 2e 


=  P11  +  P12 

where  P11  and  are  the  peak  amplitudes  of  the  two  sound  waves,  j=^l, 
^  and  are  the  propagation  vectors, 

1c  is  the  vector  from  the  origin,  and 

and  are  the  angular  frequencies  of  the  two  waves. 


1 

I 
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Then  we  have 


(pl)2  =  Pn2  +  2pllPl2  +  Pl22  . 

The  desired  sum  and  difference  frequency  components  can  be  obtained  by 
using 


<pi2>+  ♦  (Pl2). 


=  2  P11P12 


=  P  P  e 
11  12e 


+  PllP12e 


+  I  ’ll V 


where  the  +  and  -  subscript  denotes  the  sum  and  difference  frequencies. 

In  order  to  simplify  the  notation,  we  shall  use  the  subscript  (±)  on  (p^  ), 
to,  k,  and  the  absorption  coefficient  (a)  to  denote  either  the  sum  or 
difference  frequency  term.  The  upper  subscript  denotes  sum  frequency  and 
the  lower  subscript  denotes  the  difference  frequency. 

B.  Solution  of  the  Second-Order  Wave  Equation  for  Two  Plane  Waves 

We  shall  now  consider  the  Interaction  of  two  plane  waves  propa¬ 
gating  with  an  angle  9  between  their  directions  of  propagation.  For  the 
parametric  receiving  array,  we  shall  assume  that  the  first  plane  wave  is 
the  high  frequency  pump  wave  and  the  second  plane  wave  is  the  low 
frequency  signal  wave.  We  shall  ignore  the  effects  of  absorption.  The 
first-order  sound  field  is  then  given  by 


P1  "  P11  +  P12 


J  0^'  j>c&jt )  j  (1^ .  r-(n  t ) 

=  Plle  +  PI2e 


(5.1) 
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The  desired  sum  and  difference  frequency  components  of  p^  then  are 


(pl  )±  =  pnp22' 


Following  the  procedures  of  Westervelt,  *  we  find: 


16,17 


E?  (Pl2)±  to  be 

2  ,  2v  2,  2\  1  ^P1  ^ d 

IT  (Px  )±  =  V  (Px  )±  -  — g 


(o)^  )4 


c  -  dt‘ 
o 


-  I^|J 


(Pl2), 


=  ±2(k1kg-'?1*^)(p12)± 


2 

’/e  can  solve  for  (p^  )  to  get 


(A. 


2(k1k2-k1*k2) 


for  k-jkg  #  k^kg 


n 

The  time  derivative  of  (p^  )±  is 

^(pl2)±  ^  (Pi2)± 


=  T 


&  2(k1k2-¥1.^,) 


(3.2) 


(3.3) 


(3.10 


(3.5) 


! 
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and  the  wave  equation  (Equation  2.31)  is 

cfp  ,  *±S(h2)± 

^  L  2(klV1?l‘^)poc< 


(3.6) 


^1  Pl*±  1  ,  2>  1  -*  , 

— —  - ; — i  <pi  >±  -  a  “o'vVt 

c  2p  c 

ooo 


Since  can  be  expressed  as 


jCk^r-oXjt)  j(hg*r-o^t) 

P1  =  Pue  +  V 


we  have  the  following  first-order  relationships 


*i  *  -  f  /p i at  • 


(3.7) 


po^l  =  J 


P11  J^l’^V)  P12  jOyr-afet) 

— —  e  + e 

“L  “fe 


V*1 

P11^1  jO^r-cu^t)  P12%2  j(^.r-o^t) 


po®l 


(3-8) 


(3-9) 


and 


jOyr-aijt) 


P1  =  -  Ja)l  Plle 


P12e 


jO^r-a^t) 


(3-10) 


Then  the  appropriate  sum  and  difference  frequency  terms  can  be  expressed 
as 


jLtd-  \  +  ^ 

c  2  ^1P1  ±  -  2  u )iy 


P11P12  J(i4*r-a)±t) 


P  c 
0  o 


(3.11) 


1  ,  2\  "P11P12  J(k±*r-a3±t) 

—  (Pl  )±  - 


2poco 


(3.12) 


2p  c 
o  o 


ST 


-  2  P0(VU1>± - - e  ' 

By  substituting  Equations  (3.1l)>  (3*12),  and  (3*13)  into  Equation  (3*6), 
we  get 


rfp2  =lf 


2 


(pi2>± 


P  c 
Ho  o 


(3.1*0 


or 


tfpo  =  rf 


o>. 


2/r> 


(h'V 


U>±  (kjl^) 


[  20^0^ 1^)  2uxLa^(k1k2-k1«k2) 


feW 


(l  1^1 
•\2'2  *&) 


(Pl2)t 


(3-15) 


P  c 
Ho  o 


Equation  (3*15)  has  the  same  form  as  we  would  have  obtained  if  we  had  used 
Equation  (2*32)  instead  of  Equation  (2. 31)  aB  the  starting  point  for  the 
derivation.  Hence,  for  two  plane  waves  Equation  (2.31)  and  Equation  (2.32) 
are  equivalent.  The  second-order  pressure  is  then  found  to  be 


= **±{*1-% + yjfe i <Pi2>±  a _ i 

2  2a31a^(k1k2-^1.l^)  \  2  klV  PQc 


(3-16) 


This  solution  corresponds  to  the  solution  for  plane  waves  of  infinite 
extent  with  no  absorption  or  finite  amplitude  effects  present.  As  can  be 
seen,  the  second-order  sound  pressure  is  proportional  to  the  local  value 
of  (p^  )  except  when  l^k^k^kg.  This  result  led  Westervelt  to  conclude 
that  no  scattering  of  sound  by  sound  existed  when  two  sound  waves  were  not 
propagating  in  the  same  direction.  A  conclusion  based  on  Equation  (3*16) 
is  subject  to  the  conditions  that 

l)  no  absorption  is  present,  and 


28 


2)  the  plane  waves  are  of  infinite  extent. 

Strictly  speaking,  conclusions  about  the  scattering  of  sound  by  sound 
outside  the  region  of  interaction  cannot  be  made  since  the  plane  vaveB 
are  assumed  to  be  present  over  all  space. 

If  one  plane  wave  is  present  over  all  space  and  the  second 
present  only  in  the  half  space  x  g  0  and  propagates  in  the  +x  direction, 
then  we  have 


"  Plle 


j^x-o^t) 


x  g  0. 


(3.17) 


We  assume  that  the  second  wave  propagates  with  an  angle  6  with  reBpect  to 
the  direction  of  propagation  of  the  first  wave  such  that  we  have 


■12 


=  Pl2e 


jO^.r-ojjt) 


(3-18) 


If  we  neglect  the  effects  of  the  real  source  at  x  =  0,  the 
pressure  at  x  =  0  must  be  equal  to  the  local  pressure,  the  second  term  in 
Equation  0.16).  Any  solution  which  satisfies  the  second-order 
homogeneous  wave  equation 

rf $  =  0  (3.19) 


can  be  added  to  the  solution  to  satisfy  the  boundary  conditions.  One  such 
solution  is  a  plane  wave  such  that 


PllP12exp 


( kl±k2  )\ 


2(kiVVV  po°< 


where  i  is  unit  vector  in  the  x  direction, 

7  is  a  unit  vector  in  the  y  direction,  and 

kgy  is  the  component  of  the  kg  wave  vector  in  the  y  direction. 


•  r-oo±t 


(3.20) 


29 


This  solution  represents  a  plane  vave  of  frequency 

<*lVco 

2fl 


with  a  wave  vector 


We  can  express  this  solution  as 


•±  =  V 


where 


/  kp  sin  0 

F  =  -  exp  J  (kjlkg)  x  fl - —  ±  kgy  sin  0  -  r  •  (^±^) 


(5.21) 


(3.22) 


(3.23) 


(3-24) 


I  sin  0 

*  -  exp  J  (k  ±kp)  x  fl - 5-  -  k  : 

I  1  *  V 


-  kjX^kg  cos  e 


V  ' 


“±2  +  sf)  pnpia  exp  jj(^)-r-“±tj 


2(klV1V^>)  »oc< 


(5-25) 


where 


kg  =  kg  cos  0  i  +  kg  sin  0  J  ,  and 


”)  7*  “7> 

r  =  xi  +  yj  . 

If  k^  »  kg,  we  can  approximate  F  as 

/  i  kg2  sin2  @\ 

F  s  -  exp  j  (k^kg)  x  1 1  -  - - —  I  -  k^xTkg  cos  0 

L  \  / 

2  2  11 
,  kg  sin  0 

=  -  exp  i  ±kg  x  (1-cos  0)  -  g  ^  ^ — 


(3-26) 
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The  second-order  pressure  at  (x,0)  can  be  expressed  as 


p2"  =  (1+F)  p2' 


1  -  exp 


(5-27) 


.  kg2  Bin2  0 

( W  * ,  A  -  TT cos  9 


(kl±k2^ 


PllP12eXp 


2(kXk2-^L*^)  poC 


¥ - 


If  k^  »  kg,  we  have 


P2"  =  < 


1  -  exp 


k_2  sin2  0 
±  kgxd-cos  0)  -  ■ 


/ 

^1%  b\ 

*  ■ 

/_*  v 

7\  Vb  '  2A/ 

^11^12 6Xp 

j^k^xtkg  x  cos  0  ±  kgy  sin  0  -o>±tj 

■  a 

(5-28) 


2(k1k2-k1»l^  cos  0)  pQcc 


2\ _ _ _  B\__„  „  sin 


J(cn±  )(cos  0  +  2a)xP11Pj2 


2  2 
kg  x  sin  0 

V(1-=os  9)T  ur^y 


2Wo 


kgX(i-cos  0) 


exp 


'  2  2 

Lx  sin  0 

j(k1±k2)x?k2x(  1-cos  0)  -  ]^k'  —  ±  k^y  sin  0  -o>±t 


J(co±)(cos  0  +  —JPuP^x  sin [kg  ^Lk|£S_£ij 

kgX.(  1-cos  0) 


2p  c  - 
o  o 


(5.29) 


exp  j 


2  2 
kg  x  sin  0 

(VV^xd-cos  8)  -  TTi^a^r-  * 


sin  0  -0D±t 
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The  solution  in  Equation  (3.27)  gives  the  second-order  pressure  for  two 
plane  waves  interacting  over  a  half  space.  The  plane  wave  with  frequency 
f ^  propagates  in  the  x  direction  from  x  =  0. 

The  function  (1  +  F)  has  a  zero  for  0  =  0  which  cancels  the 
singularity  in  p^ ' . 

The  pressure  in  Equation  (3 .29)  has  an  amplitude  which  is  pro¬ 
portional  to  x  when  0=0.  This  result  is  consistent  with  the  result 
predicted  by  Lamb,  ^  Berktay,^0  and  Lauvstad  and  Tjotta.^  The  directivity 
function  is  the  same  as  found  by  Berktay  for  a  parametric  receiving  array 
with  a  well  collimated  pump  wave.  For  a  given  value  of  x  the  pressure 
has  a  phase  dependence  on  y  of  i  L  sin  0. 

The  second-order  pressure  is  seen  to  be  cumulative  even  when 
0*0.  Therefore,  we  see  that  cumulative  interaction  does  exist  for  an 
arbitrary  angle  of  0. 

C.  Solution  of  the  Second-Order  Wave  Equation  Using  the  Freefleld 
Green's  Functions 

Since  all  of  the  problems  of  interest  for  parametric  arrays  do 

not  have  the  simple  geometries  of  plane  waves,  it  is  necessary  to  find  a 

more  generally  applicable  solution  for  the  second-order  wave  equation.  We 

can  integrate  the  freefield  Green's  function  over  the  volume  distribution 

of  the  source.  The  freefield  Green's  function  is  a  solution  for  an 
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unbounded  medium  of  the  equation 

□  g(?Q|r)  =  -  6(ro-r)  (3. 30) 

where  r^  =  (xQ,  yQ,  zq)  is  the  observer  point, 
r  =  (x,  y,  z)  is  the  source  point,  and 
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&(?  ~"r)  =  &(x  -x)  6(y  -y)  6(z  -z)  is  the  delta  function  for  ti^ree 

O  'O  O  O 

dimensions.  Thus,  g(?o  |?)  is  the  spatial  solution  for  a  simple  harmonic 
point  source  at  r  and  has  the  form 


ikR 

8(?ol?)  =  TBT 


(3.31) 


where  R  =  I? -r  I 
'  o  1 

When  absorption  is  present,  the  solution  for  Equation  (3*30)  can  be 
modified  by  making  the  wave  number  complex.  We  then  have  a  freefield 
Green's  function  of  the  form 


(jf  -a)R 

=i4nR - 


(3.32) 


where  a  is  the  absorption  coefficient. 
A  particular  solution  for  the  equation 


can  be  expressed  in  terms  of  the  Green's  function  and  has  the  form 


(3-33) 


(3.3^) 


where  V  is  the  volume  of  integration, 

?  is  a  vector  from  the  origin  to  the  observation  point,  and 
r  is  a  vector  from  the  origin  to  the  point  of  interaction. 

D.  A  Solution  for  An  Omnidirectional  Pump  Transducer 

The  geometry  used  for  the  solution  for  an  omnidirectional  pump 
transducer  is  shown  in  Pig.  3*1*  The  signal  source  is  assumed  to  be  far 
from  the  parametric  receiver  so  that  the  low  frequency  wave  approximates 


FIGURE  3.1 

A  PARAMETRIC  RECEIVING  ARRAY 
WITH  A  POINT  SOURCE  PUMP 


34 


a  plane  wave  in  the  vicinity  of  the  receiver.  Then  at  the  point  (x,  y,  z) 
and  a  time  t  the  sound  waves  can  be  represented  in  complex  form  as 


PUe*P  [-(«!“ J^i)  \/x2+y2 +z2  -  jojjt] 

=  J - 17^ - 


(3.35) 


'+y  +z 


and 


P2  =  P12exp[-(a2-jk2)(cos  6x  +  sin  0y)-Ja^tj  (3.36) 

where  a^,  and  k^,  kg  are  the  absorption  coefficients  and  the  wave 
numbers  at  the  frequencies  ox^  and  uv,, 

P^  is  the  amplitude  of  the  pump  pressure  referred  to  1  meter,  and 
Pj2  is  the  amplitude  of  the  plane  wave  signal. 

The  direction  of  propagation  of  the  signal  wave  is  assumed  to  be 
parallel  to  the  (x,  y)  plane. 

2 

We  find  the  sum  and  difference  frequency  components  of  p^  to  be 


(Pi  )±  =  J 


/~S — ? — cT 

pnpi2exp  “ai''x  +y  +z  -ct2^x  008  0  +  y  Bin  0) 


\Z? +y2 +z2 


(3.37) 


•  exp|j(k^  >/x2+y2+z2  ±  kg  x  cos  0  ±  kgy  sin  0)  -  j(oXj±u^)tj 
The  source  term  is  then  given  by 


U+o?)  pi  i  PioexP  (-°!1  ^x2+y2 +z5  ) 


B  >2  2 

.  liHii  _A.  =  J(a,  ‘vv . 

pc  Ct2  D  c  ^  7F7Z 

00  00  * 

•  exp^-a2(x  cos  0  +  y  sin  0)J 

*  exp^j^  >/x2+y2+z2  ±  kgX  cos  0  ±  kgy  sin  0)j 

>£-  tj 


(3.38) 


exp 
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The  Integral  Green's  function  solution  for  the  second-order  wave  equation 
Is  then 


GO 


.  exp[+  >/(xo-x)^4y2+z2  -  j(cD1±o^)tj  -ixiyd* 

where  a±  Is  the  absorption  of  the  sum  or  difference  frequency  wave,  and 

xq  is  the  position  of  the  observer. 

This  rather  complicated  integral  does  not  lend  itself  to  direct 

integration.  For  its  evaluation  we  shall  use  an  alternate  method, 

commonly  used  in  optics,  which  is  known  as  the  method  of  stationary  phase. 

Single  integrals  are  frequently  evaluated  using  this  integral  technique; 

however,  double  integrals  can  also  be  evaluated  by  the  method  of  stationary 
59 

phase . 

The  solution  for  the  integral  in  Equation  (3*39)  Is  given  in 
Appendix  I.  To  the  author's  knowledge,  this  use  of  a  two-dimensional 
stationary  phase  solution  is  the  first  for  a  nonlinear  acoustics  problem. 
Thus,  the  integral  over  y  and  z  can  be  expressed  in  the  following  form 


// 


8  eJ  *  dyiz 


(3.40) 


This  integral  can  be  evaluated  by  the  method  of  stationary  phase  provided 
S(y,z)  is  a  slowly  varying  function  of  y  and  z,  and  k  f(y,z)  is  a  rapidly 
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varying  function  of  y  and  z.  The  contributions  to  the  asynrtotic  expansion 
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of  the  integral  come  only  from  certain  critical  points.  The  critical 
point  of  interest  for  us  is  the  one  within  the  region  of  integration  where 


df  _  df 
3y  “  ciz 


0 


(3-M) 


This  point  represents  the  point  on  the  wavefront  of  the  second-order 
pressure  at  the  interaction  region  which  will  eventually  strike  the 
receiver.  Fran  Appendix  A,  we  find  that 

||  =  0  f or  z  =  0  (5.42) 


and 

Tkg  sin  6  (x)(x  -x) 

^  -  0  for  y  =  k1(xQ)  ±  kgU^x) 

where  we  have  assumed  that  »  kg.  We  see  that  the  coordinate  of  the 
critical  point  varies  as  a  function  of  x. 

The  second- order  pressure  is  found  to  be 


-<W  <1+2-f>  P11P12 


p*  = 


2p  c  " 
o 


exp(jk±x 


0-jv> 


(3.^3) 
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In  the  derivation  of  Equation  (3‘^5)>  we  have  neglected  the  phase  term 


1  *22  aln2  0  /xS  \ 


(}.W 


This  term  has  a  maximum  magnitude  of 


2  2 
x  kg  sin  0 

5  Tyig"  xo  at* 


If  we  restrict  the  maximum  allowed  phase  error  to 


n 


5 


we  have 


2  2 
sin  0 

('Vkg )"  <2n 


(3.^8) 


(3^9) 


This  assumption  could  be  removed  from  the  integral  answer  by  completing  the 

2 

square  In  the  exponent  for  x  and  x  and  making  a  change  of  variables.  The 
result  would  be  in  the  form  of  a  Fresnel  integral  which  is  tabulated.  This 
solution  will  not  be  included  in  the  present  analysis  since  the  assumption 
in  Equation  (3-^9)  will  usually  be  valid  for  the  cases  of  interest  to  us. 

The  second-order  pressure  in  Equation  (3-^)  has  a  magnitude  of 


A 

P  = 


(®±)  (l+gf)  PnPi2  exp 

2p  c  5 
o  o 


-t°j ■)  *o] 


(5.50) 


and  a  directivity  function 

D(0) 


sin 

'kg  xQ(l  -  cos  0)' 

2 

K 

c 

5  Xq(1  -  cos  0) 

(3.51) 
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Except  for  absorption,  the  amplitude  is  independent  of  range, 
directivity  function  is  the  same  as  that  found  by  Berktay. 
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IV.  THE  SECOND-ORDER  SOUND  FIELD  FOR  SEVERAL  CONFIGURATIONS 
OF  THE  PARAMETRIC  RECEIVING  ARRAY 

We  shall  new  find  the  second- order  sound  pressure  for  some 
simple  pump  configurations  using  the  point- source  solution  found  in 
Chapter  III.  Since  we  are  treating  "small  signal  nonlinear"  or 
"quasilinear"  interactions,  ve  can  use  the  principle  of  superposition  on 
either  of  the  first-order  sound  fields.  In  this  case,  we  consider  the 
interactions  of  a  plane  wave  with  an  infinite  number  of  spherical  waves 
from  point  sources  summed  over  the  surface  or  volume  of  the  pump  source. 
Using  this  procedure,  we  have  a  Bimple  means  to  solve  parametric  receiving 
array  problems  for  a  wide  variety  of  geometrical  configurations.  We  shall 
only  consider  the  case  when  the  low  frequency  wave  is  a  plane  wave.  Other 
studies  by  Berktay  and  Shooter,  and  Rogers  et  al. ,  have  considered 
cases  in  which  the  wavefront  of  the  low  frequency  wave  is  spherical, 
because  the  signal  source  is  placed  a  finite  distance  from  the  parametric 
receiving  array.  In  this  analysis,  we  shall  study  only  planar  signal 
sound  waves. 

We  shall  first  consider  a  finite  length  line  source  pump 

transducer.  Two  dimensions  of  the  line  source  are  small  compared  to 

the  signal  wavelength,  and  the  other  dimension  is  on  the  order  of  a 

signal  wavelength.  The  theoretical  and  experimental  results  can  be  easily 

compared  because  the  effects  of  shadowing  of  the  low  frequency  signal  wave 

by  the  pump  transducer  are  minimized.  A  coordinate  system  which  is  the 

same  as  that  used  in  the  following  chapters  for  experimental  work  will  be 

used.  We  shall  also  investigate  the  effects  of  misalignment  of  the  line 
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source  transducer  on  the  parametric  receiving  beam  patterns.  This 
analysis  will  also  enable  us  to  investigate  the  problem  of  the  scattering 
of  sound  by  sound.  We  will  be  able  to  investigate  the  directional 
properties  of  the  second-order  sound  waves.  This  geometry  will  provide 
us  with  a  means  to  observe  the  so-called  "doppler  angles"  of  the  propaga¬ 
tion  of  the  sum  and  difference  frequency  sound  waves. 

In  Section  B,  we  shall  consider  the  use  of  a  line  hydrophone  to 
receive  the  second-order  pressure  wave  generated  by  a  point  source  and  a 
plane  wave.  The  line  hydrophone  will  have  one  dimension  which  is  on  the 
order  of  a  wavelength  at  the  signal  frequency  and  two  dimensions  which  are 
small  compared  to  the  signal  frequency  wavelength.  In  fact,  for  the 
experimental  work,  the  line  transducer  used  for  the  hydrophone  will  be  the 
same  one  which  is  used  as  a  pump  in  the  line  source  pump  experiment. 

Again,  this  geometry  will  enable  us  to  study  the  effects  of  a  misaligned 
receiver,  and  also  to  observe  the  properties  of  the  second-order  sound 
field.  We  can  observe  the  "doppler  angles"  of  the  sum  and  difference 
frequency  waves.  Again  we  shall  use  a  geometry  which  coincides  with  the 
geometry  used  for  experiments.  In  each  of  the  solutions  in  Sections  A  and 
B,  the  effects  of  the  interactions  at  the  real  source  are  not  included. 

We  simply  use  the  scattering  integral  solution  over  the  volume  of  inter¬ 
action.  We  shall  likewise  ignore  the  local  pressure  at  the  receiver. 

In  Section  C,  the  interaction  of  a  plane  wave  and  a  pump  wave 
generated  by  a  pump  transducer  which  has  all  three  dimensions,  which  are 
small  compared  to  the  signal  wavelength,  will  be  considered.  We  shall  not, 
however,  restrict  ourselves  to  the  consideration  of  pump  transducers  with 
beam  patterns  small  compared  to  the  desired  parametric  receiving  array 


beamwidth.  The  small  pump  transducer  is  of  Interest  to  us  because  such  a 
pump  transducer  vould  be  used  In  most  practical  applications.  With  a  small 
pump  transducer,  the  ratio  of  array  length  to  pump  transducer  size  is 
maximized.  The  solution  obtained  will  be  valid  for  a  range  of  pump 
transducers,  varying  from  very  narrow  beam  to  very  broad  beam.  In  each 
case,  the  observer  is  assumed  to  be  in  the  farfleld  of  the  pump  trans¬ 
ducer.  However,  the  second-order  pressure  generated  in  the  nearfield 
region  of  the  transducer  will  be  included  in  our  solution.  Likewise,  the 
solution  for  a  finite  length  line  receiver  will  be  valid  only  when  the 
omnidirectional  pump  transducer  is  in  the  farfield  of  the  receiving  trans¬ 
ducer.  A  solution  could  be  generated  for  the  near field  of  the  pump  or 
receiving  transducer  using  techniques  similar  to  those  used  by  Stenzel,^ 
Freedman,^1  and  others.^2 
A.  Solution  for  a  Line  Pump  Source 

The  solution  for  a  line  source  is  obtained  by  integration  of  the 
point  source  solution  over  the  length  of  the  line.  The  geometry  used  for 
this  solution  is  shown  in  Fig.  4.1.  The  line  of  length  2a  is  assumed  to 
lie  in  the  (x,y)  plane.  The  propagation  vector  for  the  low  frequency  wave 
is  also  in  the  (x,y)  plane. 

^  Then  the  first-order  sound  field  at  due  to  an  elemental  length 

6 1  at  ~r=£  of  the  line  source  can  be  expressed  as 

JP 

BPi (?)  =  §f  exp^ctjR-JkjR)  +  JcUjtj  (4.l) 

where  is  the  axial  farfield  pressure  referred  to  1  meter  due  to  the 
line  source.  Along  the  x  axis,  the  signal  wave  can  be  expressed  as 


FIGURE  4.1 

A  PARAMETRIC  RECEIVING  ARRAY 
WITH  A  LINE  SOURCE  PUMP 


P2  =  PI2  exPpa2“Jk2^r  COB  6  + 


(^.2) 


This  plane  wave  must  be  expressed  In  terms  of  the  new  coordinate  system 
(r,0)  which  Is  rotated  by  0"  and  translated  by  (i,0')*  The  position 
vector  in  the  new  coordinate  system  is 


(4.3) 


ir  «=  [cos  0"  (x-/  sin  0')  -  sin  0"  (y-£  cos  0*)]1* 

+[sin  0"  (x-i  sin  0')  +  cos  0"  (y-£  cos  0')!^ 

where  1*  and  J  cure  the  basis  vectors  for  the  new  coordinate  system.  The 
plane  wave  vector  in  the  new  coordinate  system  is 

cos (0"+0)J  +  kg  sin(0"+0)]£  (4.4) 

and  then  T^,*r;  is  given  by 

*  r  -  kg/  sin(04€')  .  (4.5) 

The  plane  wave  pressure  expressed  in  terms  of  r  ,9,9' ,  and  0"  is 

P2  *  cos (0*0 ")  +  i  sin(0+0')]}  .  (4.6) 

Then  from  Equation  (3-46)  the  sum  and  difference  frequency  sound  pressure 
at  L  is 


P  P 
*  ^  ^ 


5p(L,0,0',0")  =  -  f  12 

2a(2p  c n 


o  o 


&/  exp|-ag/  sin(0+0*)  -  (a±-Jk±)R| 

(4.7) 


•  expjijkg/ 


sin(0+0 ' )  T  j  —5—  (l  -  cos(0+0 ") ) 


kgR 

~2 


sin!-~-  (l  -  cos(0+0")) 


kpR 

-f-  (1  -  cos(0+9")) 


■n 

where  p  =  1  +  t-t-  . 
r  2A 
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Using  the  farfield  approximation  for  a  line  source,  we  have 

„  L-i  sin  S'  ^ , 

cos  6 

for  the  phase  and  amplitude  terms  of  the  elemental  pressure.  Ihen 
elemental  pressure  can  be  shown  to  be 


Bp(L,0,0')  = 


2a(2poco3) 

exp[-(a±-Jk±)l]exp[jJ>^JAi-Ja>±tJ  Bl^iC> 


where 


kp 

A  =  ±  ~  sin(0*0')  ±  gS.  sin  0*  -  k±  sin  0' 

kp 

C  =  [-sin  0’  +  sin(0+0’)]  ,  and 

kp 

M  =  ^  [L(l  -  cos  0)] 


The  total  pressure  can  be  obtained  by  integrating  over  i  such  that 

p(L,0,0’)  =  ^  exp[-(a±-jk±)L]  .  ^  J 


2poco 


exp[*JM  +  JAi]  di 

If  we  assume  that  sin(M<-iC)/m-iC  is  a  slowly  varying  function  of  i, 

,,  „  "^a)l±afe^PllPI2  T  /..  .  .1  sin  aA  Bin  M 

p(L,0,0  )  s - 5 - exp|-(a±-jk±)I*jM-Jo>±tJ  — ^ - — 


2p  c 
o  o 


(4.8) 
the 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 


.  (4.14) 
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The  maximum  amplitude  for  6=0  Is 


P  "  3 

2p  c  ; 
K0  0 


exp[-(a±L)] 


(4.15) 


The  approximation  used  in  Equation  (4.l4)  can  be  eliminated  by  direct 
integration  of  Equation  (4.13). 

The  pressure  can  be  expressed  as 

P  =  exp £+jk±I^ JM- jco±t j  •  J  exptjAi]  .  (4.l6) 

[exp[j(MfCf)]  -  exp[-j(M+Ci)]l 

2j(M+CiJ  * 


and  vith  a  change  of  variables,  we  have 


p  ■  Tsfer  1  P  *  -  P  * 

yl  y3 

=  i2sfey  exp[+Jk±L?JM  -  -4Ci(y2)  ■  Ci(yi)  -  Ci(y4)  +  Ci(y3)| 

+  [si(y2)  -  Si(y]L)  -  Si(y4)  +  Si(y3)j  (4.17) 

where  Ci(y)  is  the  cosine  integral,  Si(y)  is  the  sine  integral, 


yl  =  (M"aC) 

y2  =  ^£l  (M+aC) 

y3  =  (M_aC) 

y4  =  ^lP'  (M+aC^ 


9 


and 


(4.18) 
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In  order  to  evaluate  the  cosine  Integral  In  the  case  vhen  a  pair 
of  the  yi  's  are  zero,  the  limit  as  each  approaches  zero  must  be  taken 
in  order  to  obtain  a  solution.  Likewise,  if  any  y  equals  0  over  the  range 
of  Integration,  Cauchy's  integral  formula  must  be  used  to  obtain  the  value 
of  the  integral. 

Computed  results  using  Equation  (4.17)  will  be  presented  in 
Chapter  VI.  The  basic  properties  of  a  parametric  receiving  array  using 
a  line  source  pump  can  best  be  seen  by  an  examination  of  Equation  (4.l4). 
The  amplitude  is  seen  to  be  independent  of  range  except  for  absorption. 
This  result  is  subject  to  the  condition  that  we  are  in  the  farfield  of  the 
pump  transducer.  However,  the  interaction  from  the  nearfield  region  of 
the  transducer  is  included  in  the  result. 

Further  examination  of  the  expression  for  the  second-order 
pressure  indicates  that  we  have  a  directivity  function  which  is  a  product 

V 

of  two  functions;  one  function  is  ^  ,  where  M  =  L(l  -  cos  6), 

which  is  the  directivity  function  of  an  end-fire  array  of  length  L.  This 

function  is  the  same  as  the  directivity  function  for  two  plane  waves 

sin  aA 

interacting  over  some  distance  L.  The  second  function  . ,  where 

Jr  Jp 

A  =  ±  sin(0+0')  ±  g—  sin  0'  -  k±sin  0'  , 

is  an  interesting  one  because  this  function  includes  both  0  and  0'.  For 
the  case  when  0=0  the  dependence  of  the  directivity  function  on  0'  is 
the  same  as  that  of  the  original  pump  wave.  This  result  does  not  agree 
with  the  results  of  a  previous  investigation  by  Rogers  et  al.,  in  which 
the  second-order  sound  field  was  predicted  to  have  the  properties  of  the 
pump  wave  as  if  it  were  operated  at  the  sum  or  difference  frequency.  We 


48 


shall  examine  this  result  further  in  Chapter  VI  vhen  It  Is  compared  to 
experimental  data.  With  further  examination  of  the  expression  for  A  ve 
see  that  the  array  directivity  function  as  the  angle  0  is  changed  will  be 
modified  for  various  values  of  9'.  This  result  is  also  in  conflict  with 
the  result  predicted  by  Rogers  et  al.  This  result  will  also  be  compared 
extensively  with  experimental  data  since  it  allows  us  to  examine  the 
problem  of  scattering  of  sound  by  sound  when  the  two  waves  are  propagating 
at  nonzero  angles. 

When  0'  is  zero,  ve  observe  that  the  low  frequency  directivity 
of  the  parametric  receiving  array  will  be  modified  slightly  by  the 
expression 

sin  aA 
aA 


This  modification  is  less  than  one  might  intuitively  expect.  The  reduction 
in  the  side  lobes  is  less  than  that  produced  by  taking  the  product  of  the 
directivity  function  of  the  end-fire  array  and  directivity  function  of  the 
pump  at  the  signal  frequency. 

B.  Solution  for  a  Line  Hydrophone  and  a  Point  Source  Pump 

As  in  the  previous  example,  we  shall  use  the  point  source 
solution  as  a  starting  point.  In  the  solution  for  the  point  source  pump, 
a  small  point  receiver  is  assumed.  Using  the  geometry  shown  in  Fig.  4.2, 
we  shall  consider  response  of  a  line  hydrophone  of  length  2a  consisting 
of  an  infinite  number  of  point  receivers  lying  in  the  (x,y)  plane.  The 
propagation  vector  for  the  low  frequency  wave  is  also  in  the  (x,y)  plane. 

The  receiving  sensitivity  due  to  an  elemental  length  b£  at  "r  is 

M  (o>)BJ 


8M(>)=  2a 


(^.19) 


2a 
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where  M  (co)  is  the  receiving  response  of  the  line  in  volts  per  micropascal 
c 

to  a  plane  wave  of  frequency  os  propagating  in  a  direction  perpendicular  to 
the  line. 

The  high  frequency  pump  wave  has  the  form 


P11  "  J  "IT 


(4.20) 


where  is  the  farfield  pressure  referred  to  1  meter.  The  signal  wave 
at  (r,0")  can  be  expressed  as 


P12  =  P12  cos(0-e")-jo^t] 


(4.21) 


The  second-order  pressure  at  (r,0M)  is 


P2(R,e")  = 


(oijfcl^) 


2p  c 
o  o 


5  @P11PJ2 


•  expj-(a±-Jk±)  H  ±  JkgR  ii- 

«»[yi  (* 

'  y  (x  HMEEZf 


(4.22) 


Again  using  the  farfield  approximation,  ve  have 

a  .  L  '  4  ?'  s  l  -  i  sin  e' 

cos  0' 


(^.25) 


for  the  phase  and  amplitude  term.  The  elemental  voltage  generated  by  the 
second-order  sound  pressure  is 


BE  =  - 


McBi  (o3l±o^)3P11P12 


2a 


(4.24) 


2p  c 
o  o 


•  exp|-(a±-jk±)L-jajtJ 

.  expfTj  ^  [L(l  -  cos  0)]  +  JBi)  •  S1^ 
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where 


B  =  ±  ^  sin(e-e')  ±  -  k±  sin  9' 

D  =  l-sin  9'+sin(0-0')]  >  and 


M  =  ~  [L(l  -  cos  6)] 


The  total  voltage  can  be  obtained  by  integrating  over  l,  yielding 


E  - - exp[-(a±-jk±)L-jcD±t 


2p  c 
o  o 


P  P 
U  ’ 

3  c 


•  s  /a  u 

-a 


(4.25) 

(4.26) 

(4.27) 


(4.28) 


If  we  assume  that  is  a  slowly  varying  function  of  i,  we  have 

(u^±ol)PP  P  M  .  , 

E  S - L-Z - -A  ■  ■  exp[-(a±-jk±)L¥jM-ja>±tJ 


2p  c 
o  o 


sin  aB  sin  M 
*  aB  M 

The  maximum  amplitude  for  0  =  0  is 

HV^iiVc 


E 


max 


exp  [-a+L] 


pc 

Ko  o 


(4.29) 


(4.30) 


The  approximation  used  In  Equation  (4.29)  can  be  eliminated  by 
direct  integration  of  Equation  (4.28). 
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The  received  voltage  can  be  expressed  as 
E 

E  =  exp  |+jk±LrJM-jo)±t  ]  (4. 31) 

-a, 

and  with  a  change  of  variables  we  have 
E 

E  =  2a^2D)  exp[+<3k±LT«3MfJ  jp  "  Ju>±t]  (4.32) 

-  j[Ci(y2)  -  Ci(y1)  -  Ci(y4)  +  Ci(y3)) 

+  [Si(y2)  -  Si(yi)  -  Si(y4)  +  Si(y3)] 
where  Ci(y)  is  the  cosine  integral,  Si(y)  is  the  sine  integral, 

?!  *  ^  <“-»)  - 

y2  "  (M+aD) 

(4.33) 

y3  =  ‘^ET^’  >  and 

y4  =  '^TI)'‘)  (MfaD) 

The  solution  for  a  line  receiver  is  very  similar  to  the  solution 
for  the  line  pump.  However,  upon  examination  of  Equation  (4.29),  we  note 
that  there  is  an  important  difference.  The  beam  pattern  obtained  by 
varying  6'  with  0  =  0  is  the  same  as  that  obtained  by  the  line  receiver  at 
the  sum  or  difference  frequency.  When  0 '  is  equal  to  0,  the  parametric 
receiving  array  beam  pattern  is  the  same  as  that  found  for  a  line  source 
pump.  These  properties  will  be  demonstrated  in  Chapter  VI  when  we 
compare  the  experimental  and  theoretical  results. 
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C.  Kolutlons  for  the  Parametric  Receiving  Array  for  Pump  Source  with 
T\ro  Finite  Dimensions 

We  shall  now  consider  the  interaction  of  a  pump  wave  produced  hy 
a  rectangular  piston  with  one  dimension  which  is  small  compered  to  signal 
frequency  wavelength,  and  a  second  dimension  which  can  be  a  wavelength 
or  more  at  the  signal  frequency.  The  planar  radiator  with  dimensions  2a 
and  2b  is  assumed  to  have  an  angle  0*  with  respect  to  the  plane  x  =  0 
as  shown  in  Fig.  4.3.  The  observer  is  in  the  plane  z  =  0.  From 
Equation  (4.17),  the  solution  for  an  elemental  area  of  length  2a  and 
height  6z  for  z  =  0  is 

6p  =  (2b rlf  2c  |exp["(°1“*Jk±)  l*jm  -  -  jv]  | 

•  |-j[Ci(y2)  -  Ci(y1)  -  Ci(y^)  +  Ci^)]  (4.34) 

+  [Si(y2)  -  Si(yi)  -  Si(y^)  +  Si(y5)]  | 

where  y^^  is  as  defined  in  Equation  (4.l8). 

If  2b  is  small  compared  to  wavelength,  the  contributions  from 
the  various  elemental  areas  2a  6M  will  be  equal  for  a  given  6M  in  the 
farfield  of  the  radiator  since  R '  s  L  and  0 ' 1 '  s  0  for  small  values  of  M. 

The  total  pressure  is  the  SFjne  as  that  found  in  Equation  (4.17). 
If  2a  is  small  compared  to  the  signal  frequency  wavelength, 

(cn1±o^)pP11P;i2  sin  aA  sin  M 

P  ■ - j - exp|-(a±-jk±)  D:jM-ja)±t]  •  — ^ - jj—  •  (4.35) 

2poco 

We  shall  now  consider  a  circular  piston  pump  transducer  with  a 
radius  which  is  small  compared  to  a  signal  frequency  wavelength.  The 
geometry  for  this  solution  is  shown  in  Fig.  4.4.  Again  we  shall  use  the 
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solution  for  a  line  source  as  a  starting  point.  In  the  farfield  R'  is 

approximately  equal  to  L,  and  the  pressure  contribution  from  a  line  source 

n — 5" 

of  length  va  -z  and  width  of  5z  is  given  by 


8p  =  - 


2(,o°c3 


(4.36) 


•  exp 


(a±-jk±LTjM-joj±t) 


sin  M 
M 


sin 


*4 


T’a 


6z 


where 


fS2) 


is  the  normalization  for  the  area  found  by  integrating 


the  incremental  area  over  61  such  that  we  have 


(M7) 


We  now  consider  the  integral 


2  sin  'JeF-z  A 


dz 


A  II  a 


With  a  change  of  variables,  this  integral  is 


t  sin  (t  a  A)  dt 
allA  A  -  td 


Thic  Integral  can  be  evaluated  using  Sonine's  integral  described  by 
Watson.^  This  integral  is  a  definition  of  a  Bessel  function  of  order  m 
and  has  the  form 


P  m-n  A 

I  (x)  =  - — -  /  J 

ffi  2m-nr(M-n)  J0  n 


(xt)  tn+1  (l-t2)m-n-1  dt 


(4.38) 
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where  J  (x)  is  a  Bessel  function  of  order  M, 
m 

T  is  the  Gamma  function,  and 

Jn(xt)  is  a  Bessel  function  of  order  n. 

For  n  =  ^  and  M  =  1,  this  integral  has  the  form 


Using  Equation  (4.39)#  we  have 


(4.39) 


4 

(aA)ll 


t  sin(taA) 
v/l-t2 


dt 


2  Jx(aA) 
aA 


(4.4o) 


The  second-order  pressure  for  a  circular  pump  source  is  then  found  to  be 


(co  ±cn  )pp^  P  - 

l  - - exp  U(a+-jk±L)  *  jm  -  Jaj.t 

2p  c  ;  1 

*0  o 

•  w  2  J.(aA) 
sin  M  lx 

M  aA 


(4.41) 


The  amplitude  of  the  pressure  expressed  in  Equation  (4.4l)  is 
the  same  as  that  found  for  a  rectangular  piston  transducer.  The  direc¬ 
tivity  function  is  a  product  of  the  end-fire  array  function  and  the 
2  Jp(aA)  2  J  (aA) 

function  - gy  ■  .  The  function  - ^ - resembles  the  directivity 

function  for  a  circular  piston  found  in  linear  theory.  For  the  second- 
order  sound  pressure,  A  replaces  the  expression  k^  sin  0'.  The  properties 
of  the  parametric  array  receiving  with  a  circular  piston  pump  for  various 
0  and  0'  are  Been  to  be  very  similar  to  these  properties  of  a  parametric 
receiving  array  with  a  square  piston  pump  if  we  take  into  account  the 
difference  in  the  shape  of  each  aperture. 


V.  EXPERIMENTAL  APPARATUS 


The  experimental  study  of  the  parametric  receiving  array  was 
conducted  aboard  the  STEP  Barge  at  the  Applied  Research  Laboratories  Lake 
Travis  Test  Station.  This  barge,  designed  for  testing  large  sonar 
transducers,  has  heavy  handling  equipment  available  for  moving  and  rotating 
large  transducer  arrays  and  electronic  instrumentation  to  do  acoustical 
and  electrical  performance  testing  on  these  transducers. 

The  STEP  Barge  is  moored  in  100  ft  of  water  on  Lake  Travis,  a 
fresh  water  lake  about  25  miles  west  of  Austin,  Texas. 

The  mechanical  and  electrical  apparatus  used  for  the  parametric 
receiving  array  experiments  is  described  in  the  following  two  sections. 

A.  Mechanical  Apparatus 

An  overall  view  of  the  parametric  receiving  array  is  shown  in 
Fig.  5.1.  The  parametric  receiving  array  was  suspended  below  the  STEP 
Barge  using  an  aluminum  I-beam  attached  to  a  column  which  could  be  rotated. 
The  low  frequency  source  used  to  generate  the  signal  wave  was  located  on 
a  second  barge  some  525  ft  away  from  the  center  of  the  parametric 
receiving  array. 

In  Fig.  5.2,  we  see  a  close-up  of  the  parametric  receiving 
array.  The  beam  used  for  supporting  the  two  transducers  was  constructed 
so  that  the  two  ends  could  be  lifted  for  removal  or  replacement  of  the 
transducers  with  the  column  in  place  in  the  water.  Support  brackets  were 
mounted  on  either  end  of  the  I-beam  for  the  placement  of  the  two  trans¬ 
ducers.  The  separation  distance  between  the  two  transducers  was 
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FIGURE  5.1 

NONLINEAR  PARAMETRIC  RECEIVING  ARRAY  EXPERIMENT 


DETAILED  VIEW  OF  THE  PARAMETRIC  RECEIVING  ARRAY 
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adjustable.  For  the  experiments  described  In  the  next  chapter,  the 
transducers  were  placed  a  fixed  48  ft  distance  apart  and  equidistant  from 
the  support  column.  The  support  on  one  end  of  the  beam  had  a  small 
rotator.  The  angle  of  rotation  of  the  transducer  mounted  on  this  end 
corresponded  to  6'  in  the  theoretical  model.  The  I-beam  was  rotated  by 
the  column  mechanically  coupled  to  a  synchro.  The  angle  measured  by  this 
synchro  corresponded  to  the  angle  6  in  the  theoretical  model.  The  two 
transducers  for  the  parametric  receiving  array  were  mounted  6  ft  below 
the  I-beam  in  order  to  minimize  the  effects  of  acoustical  reflections. 

The  I-beam  was  supported  27  ft  below  the  water  surface,  so  that  the 
parametric  receiving  array  was  33  ft  below  the  water  surface. 

B.  The  Measurement  System  for  the  Parametric  Receiving  Array 

An  overall  block  diagram  for  the  system  used  for  the  parametric 
receiving  array  experiments  is  shown  in  Fig.  5-3*  The  pump  transducer 
was  driven  by  a  power  amplifier  with  a  CW  signal  generated  by  a  precision 
crystal  oscillator.  The  pump  frequency  remained  fixed  at  90  kHz  throughout 
all  of  the  experiments  conducted  in  this  study.  Special  care  was  taken 
to  ensure  that  the  side  band  noise  level  generated  by  the  power  amplifier 
was  low  in  order  to  achieve  a  good  signal-to-noise  ratio  at  the  side  band 
frequencies.  Figure  J.k  gives  a  breakdown  of  the  receiving  subsystem  for 
the  parametric  receiving  array.  The  signal  from  the  hydrophone  was  fed 
directly  into  a  crystal  notch  filter  providing  an  80  dB  rejection  of  the 
carrier  signal.  The  signal  was  then  amplified  and  fed  into  a  second 
crystal  filter  to  further  remove  the  carrier  frequency.  Two  passive 
bandpass  filters  were  used  to  select  either  of  the  side  bands  at  the  sum 
or  difference  frequency.  A  sampling  digital  voltmeter  measured  the  rms 


62 


I 

* 


FIGURE  5.3 

PARAMETRIC  RECEIVING  ARRAY  BLOCK  DIAGRAM 
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FIGURE  5.4 

RECEIVING  SUBSYSTEM  FOR  THE  PARAMETRIC  RECEIVING  ARRAY 
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value  of  the  side  band  signal  at  a  specified  time.  For  the  calibration 
of  the  parametric  receiving  array,  the  voltage  was  read  directly  from  the 
digital  readout.  Beam  patterns  were  plotted  on  either  a  polar  or  rectan¬ 
gular  plotter  using  the  dc  output  of  the  sampling  digital  voltmeter.  The 
dc  output  was  also  used  with  a  portable  voltohmmeter  for  the  alignment  of 
the  various  transducers  during  the  experiments. 

A  block  diagram  of  the  low  frequency  source  subsystem  is  shown 

in  Fig.  5*5*  A  gated  sine  wave  at  the  signal  frequency  was  used  to  drive 

the  farfield  transducers.  Special  low  pass  filters  were  used  on  the 

input  and  output  of  the  power  amplifier,  thereby  eliminating  noise  at  the 

side  band  frequencies  of  the  parametric  receiving  array.  The  calibration 

of  the  parametric  receiving  array  was  achieved  in  several  steps.  First, 

the  transducers  used  for  absolute  sound  pressure  level  measurements  were 

calibrated  using  the  3-transducer  reciprocity  calibration  procedure 

6k 

described  by  Bobber.  A  block  diagram  for  the  system  used  for  this 
calibration  is  shown  in  Fig.  5-6.  The  calibration  of  the  transducers 
was  made  with  the  transducers  in  place  in  order  to  minimize  the  errors 
caused  by  transducer  misalignment.  Two  additional  transducers  were  used 
to  assist  in  the  calibration:  one  was  mounted  on  a  movable  trolley  and 
the  second  was  suspended  by  its  cable  at  the  center  of  the  array.  Since 
the  crystal  filter  was  coupled  directly  to  the  hydrophone,  a  second 
calibration  was  necessary  to  correct  for  the  loading  by  the  crystal 
filter.  This  calibration  was  accomplished  using  the  comparison  method. 
First,  the  open  circuit  voltage  of  the  hydrophone  was  measured  for  a 
given  acoustic  pressure;  then,  the  output  voltage  of  the  hydrophone,  filter, 
and  amplifier  combination  was  measured  at  the  various  side  band  frequencies. 
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The  low  frequency  signal  pressure  amplitude  was  measured  using 
a  low  frequency  standard  hydrophone,  also  calibrated  using  the  reciprocity 
calibration  procedure.  The  low  frequency  signal  had  to  be  measured  with 
great  care  because  the  signal  propagated  over  a  long  distance  to  the  para¬ 
metric  receiving  array,  and  had  to  be  measured  early  in  the  pulse  before 
reflections  from  the  surface  arrived.  On  the  other  hand,  the  mechanical 
Q  of  the  source  transducer  increased  the  rise  time  of  the  signal  to 
several  periods.  Consequently,  the  number  of  cycles  available  for  making 
a  measurement  was  small.  To  ensure  a  good  measurement  of  the  low 
frequency  signal,  the  sampling  digital  voltmeter  was  used  in  the  FER 
INTEGRAL  RMS  mode.  This  mode  makes  the  sample  time  an  Integral  number  of 
periods  at  the  signal  frequency.  The  timing  for  this  measurement  is 
derived  from  a  reference  signal  from  the  low  frequency  oscillator.  The 
low  frequency  standard  hydrophone  was  placed  immediately  in  front  of  the 
parametric  receiving  array  as  is  illustrated  in  the  overall  system 
diagram  (Fig.  5.1). 

While  great  care  was  used  to  make  the  acoustical  measurements, 
it  should  be  mentioned  that  errors  as  great  as  12  dB  could  occur  on  any 
particular  measurement  because  of  uncertainties  in  the  conditions  in  the 
medium.  In  our  experiments,  an  attempt  to  improve  this  accuracy  to  better 
than  12  dB  was  made  by  calibrating  the  various  transducers  in  situ  and 
by  measuring  the  rms  value  of  the  signals  in  integral  number  of  cycles 
using  a  sampling  digital  voltmeter  with  a  measurement  resolution  of 
±0.1  dB.  These  precautions  still  do  not  preclude  the  possibility  of 
irregularities  in  the  propagation  medium. 
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Five  types  of  transducers  were  used  for  the  parametric  receiving 
array  experiments.  They  are  as  follows. 

1.  A  rectangular  transducer  l8  in.  long  and  3  in.  high  for  a 
pump  and  a  receiver  transducer. 

2.  A  small  standard  transducer  with  an  omnidirectional  trans¬ 
ducer  beam  pattern  for  a  pump  and  receiver  transducer. 

3.  A  low  frequency  standard  transducer. 

4.  A  4  in.  piston  transducer  for  a  pump  end  a  receiver 
transducer . 

5.  A  broadband  standard  hydrophone. 

An  additional  description  of  these  transducers  is  given  in 


Appendix  B. 


VI.  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 
FOR  THE  PARAMETRIC  RECEIVING  ARRAY 

We  shall  now  examine  the  theoretloal  results  of  Chapters  III  and 
IV  In  light  of  some  experiments  on  the  parametric  receiving  array.  The 
parametric  receiving  array  for  these  experiments  had  a  fixed  length  of 
48  ftj,  a  limitation  which  was  minimized  by  using  various  signal  frequencies 
and  various  transducers  for  our  experiments. 

Unfortunately,  one  cannot  measure  the  interaction  of  two 
infinite  plane  waves  experimentally  to  compare  to  the  theoretical  model 
of  Chapter  III.  We  can,  however,  observe  experimentally  a  spherically 
spreading  wave  generated  by  a  point  source  interacting  with  a  plane  wave. 
Perhaps  this  case  represents  the  physical  situation  in  which  the  theory 
might  best  predict  the  experimental  results.  We  shall  likewise  examine 
the  performance  of  the  parametric  array  with  a  rectangular  piston  pump 
with  one  dimension  on  the  order  of  a  wavelength  at  the  signal  frequency 
and  one  dimension  that  was  small  compared  to  a  wavelength  at  the  signal 
frequency.  This  case  is  interesting  because  it  allows  us  to  look  at  the 
preferred  angle  of  propagation  of  the  second-order  sound  field  as 
predicted  by  1  theory.  The  properties  of  a  parametric  receiving  array 
using  a  line  receiver  were  measured.  A  small  piston  pump  was  also  examined 
to  compare  to  the  theory  presented  in  Chapter  TV.  In  each  of  these  cases, 
the  geometries  of  the  theory  and  experiment  were  made  to  coincide, 
thereby  facilitating  comparison  between  the  theory  and  experiment.  These 
geometries  are  shown  in  Fig.  6.1. 
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(o  )  A  PARAMETRIC  RECEIVING  ARRAY 
WITH  A  LINE  SOURCE  PUMP 


(b)  A  PARAMETRIC  RECEIVING  ARRAY 
WITH  A  LINE  RECEIVER 


FIGURE  6.1 

GEOMETRIES  FOR  THE  PARAMETRIC 
RECEIVING  ARRAY  EXPERIMENTS 
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A.  Experiments  Using  a  Point  Source  Pump 

A  small  transducer  vhich  was  omnidirectional  in  the  horizontal 
plane  was  mounted  rigidly  6  ft  below  one  end  of  the  parametric  receiving 
array  I-beam.  On  the  opposing  end  of  the  I-beam,  an  assortment  of 
transducers  were  mounted.  These  transducers  included  an  omnidirectional 
transducer,  a  small  piston  transducer,  and  the  rectangular  transducer. 
Parametric  receiving  array  beam  patterns  at  5  kHz  for  the  omnidirectional 
pump  transducer  are  shown  in  Fig.  6.2.  These  beam  patterns  were  made 
using  the  4  in.  square  piston  transducer  as  the  receiver.  Additional 
beam  patterns  using  the  line  transducer  as  a  receiver,  with  the  omni¬ 
directional  pump  transducer,  will  be  shown  in  Section  VI. C.  since  our 
theoretical  result  showed  that  a  large  receiving  transducer  would  modify 
the  beam  pattern. 

The  result  from  Equation  (3-46)  was  uBed  to  calculate  the 
theoretical  beam  pattern,  and  an  HP  9^50  computer  program  was  written  to 
plot  this  result.  This  program,  listed  in  Appendix  C,  was  also  used  to 
plot  the  theoretical  beam  patterns  for  the  rectangular  pump  transducers. 
As  can  be  seen,  the  theoretical  model  does  predict  the  correct  beam 
pattern  for  the  parametric  receiving  array.  The  major  lobe  width  and 
side  lobe  height  and  bearing  were  very  close  to  that  predicted  by  the 
theoretical  model.  This  remarkable  agreement  between  the  theory  and 
experiment  leaves  little  doubt  that  this  model  is  accurate  for  predicting 
the  performance  with  an  omnidirectional  pump.  However,  the  expression 
for  the  phase  of  the  omnidirectional  pump  would  not  necessarily  be 
verified  by  the  excellent  agreement  of  the  amplitude.  Amplitude  measure¬ 
ments  were  made  for  the  omnidirectional  transducer  using  the  procedures 


(a)  SUM  FREQUENCY 


(b)  DIFFERENCE  FREQUENCY 


FIGURE  6.2 

PARAMETRIC  RECEIVING  ARRAY  BEAM  PATTERNS 
USING  AN  OMNIDIRECTIONAL  PUMP  AT  5  kHz 
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described  in  Chapter  V.  Experimental  measurements  were  also  made  using 

an  omnidirectional  receiver.  The  3  dB  beamwidth  was  found  to  be  equal  to 

that  of  the  theoretical  result.  However,  the  signal-to-noise  ratio  was 

too  low  to  measure  the  side  lobe  levels. 

The  sum  frequency  sound  pressure  level  as  a  function  of 

frequency  is  shown  in  Fig.  6.3.  The  sound  pressure  level  has  been 

normalized  to  the  level  l'or  sound  pressure  levels  of  0  dB  re  1  uPa  at 

1  m  for  the  pump  pressure  and  0  dB  re  1  pPa  for  the  signal  pressure  at 

the  receiver.  The  theoretical  curve  for  the  sound  pressure  level  was 

plotted  using  Equation  (3. ^+6)  for  6=0. 

The  amplitude  of  the  measured  sound  pressure  was  corrected  to 

take  into  account  the  finite  distance  to  the  position  of  the  low  frequency 
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source.  This  correction  was  predicted  by  Berktay  and  Shooter.  Taking 
into  account  the  spherical  nature  of  the  low  frequency  wavefront,  the 
second-order  sound  pressure  for  0=0  is 
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where  R  is  the  range  of  the  low  frequency  source  from  the  receiving  trans¬ 
ducer.  For  our  experiment,  the  range  correction  was  O.63  dB.  Berktay  and 
Shooter  also  demonstrated  that  the  discrepancies  in  beam  patterns  would 

be  negligible  if  —  >  b  and  the  array  was  rotated  about  its  center.  In  our 
L 

p 

case,  ^  was  approximately  7*  For  the  computation  of  the  sound  pressure 
level,  we  used  a  value  of  3.5  for  (l  +  ~)  as  predicted  by  Beyer, 

Coppens  et  al.,^  for  a  water  temperature  of  70°F.  Unfort unately,  as  was 
mentioned  in  Chapter  V,  the  basic  experimental  error  could  be  as  large 
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as  2  dB,  were  the  number  of  pressure  measurements  made  taken  into 
account.  However,  since  careful  measurements  were  made  at  several 
frequencies,  the  overall  accuracy  could  be  improve!  somewhat.  The 
frequency  response  of  the  sum  pressure  is  seen  tc  have  a  frequency 
dependence  of  (ca^+a^) .  Some  deviation  in  this  dependence  at  the  higher 
frequencies  occurred  due  to  a  lower  signal  to  noise  ratio  for  the 
measurement  signal  caused  by  characteristics  of  the  low  frequency  signal 
source. 

The  difference  frequency  sound  pressure  level  is  shown  in 
Fig.  6.4.  Again,  this  normalized  sound  pressure  level  haB  been  corrected 
for  the  finite  range  of  the  low  frequency  sound. 

B.  Experiments  on  a  Parametric  Receiving  Array  with  a  Line  Source  Pump 

In  the  experiment  described  now,  a  rectangular  source  was  used 
as  an  approximation  to  a  line  source  to  study  the  properties  of  the  second- 
order  sound  field.  This  rectangular  transducer  was  constructed  with 
bizonal  shading  to  reduce  the  side  lobeB.  The  velocity  of  this  type  of 
acoustical  source  is  not  uniform  across  the  aperture.  The  center  portion 
of  the  transducer  moves  with  twice  the  velocity  of  the  ends  of  the  aperture 
for  a  given  drive  voltage.  This  discrepancy  between  the  theoretical  model 
and  experiment  configuration  will  be  eliminated  by  modifying  the  theory  in 
Chapter  IV  to  include  the  effects  of  bizonal  shading.  The  sum  of  the 
results  of  two  rectangular  apertures  will  be  summed:  one  is  the  full 
width  of  the  transducer,  generating  a  given  far field  sound  pressure  level 
referenced  to  1  m,  and  the  second  transducer  aperture  is  smaller  and  has 
the  appropriate  length  and  amplitude  weighting  to  give  the  approximate 
side  lobe  reduction  of  the  active  beam  pattern.  The  amplitude  of  the  sum 
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of  the  two  sound  pressures  was  adjusted  to  be  P1  referenced  to  1  in,  the 
value  used  In  the  theoretical  result  (Equations  4.l4  and  4.17).  The 
velocity  amplitude  of  the  center  portion  of  the  aperture  was  adjusted  so 
that  it  had  twice  the  velocity  of  the  outer  portion  of  the  aperture.  The 
fact  that  the  correction  for  the  shading  cap  be  easily  made  demonstrates 
the  versatility  of  the  theory  that  has  been  developed,  dfce  beam  pattern 
for  the  narrowbeam  pump  transducer  at  90  kHz  is  shown  in  Pig.  6.5.  The 
results  are  plotted  in  rectangular  coordinates  to  enhance  the  comparison 
between  the  theory  and  the  experiment.  The  beam  pattern  was  obtained  by 
varying  0'  and  observing  the  pump  signal.  The  bizonal  shading  of  the 
transducer  causes  the  usual  beam  pattern  to  be  somewhat  modified  with  a 
broader  beam  and  lower  side  lobes  for  a  given  aperture.  In  this 
particular  example,  an  outer  aperture  of  I.38  ft  and  an  inner  aperture  of 
O.78  ft  was  used. 

An  oscilloscope  waveform  of  the  sum  and  difference  frequency 
side  bands  observe  \  at  the  output  of  the  second  crystal  filter  is  shown 
in  Fig.  6.6(a).  The  signal  frequency  is  5  kHz  and  the  pump  frequency  is 
90  kHz.  The  pulse  envelope  of  the  difference  frequency  side  band  signal 
is  shown  in  Fig.  6.6(b). 

The  sum  and  difference  frequency  sound  pressure  levels  as  a 
function  of  frequency  are  shewn  in  Fig.  6.7  and  Fig.  6.8.  Again,  the 
sound  pressure  level  has  been  normalized  for  the  sound  pressure  level  of 
0  dB  re  1  pPa  at  1  m  for  the  pump  pressure  and  0  dB  re  1  (iFa  for  the 
signal  pressure  at  the  receiver.  The  correction  for  the  finite  range  of 
the  low  frequency  source  was  again  incorporated  in  the  result.  The 
agreement  between  the  theory  and  experiment  was  excellent.  A  value  of  5 
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FIGURE  6.5 

LINE  SOURCE  PUMP  BEAM  PATTERN 
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was  used  for  B/A.  Difference  frequency  begun  patterns  of  the  parametric 
receiving  girray  at  signal  frequencies  of  3  and  5  kHz  are  shown  in  Fig.  6.9. 
Sum  frequency  begun  patterns  were  found  to  be  very  similar  to  the  difference 
frequency  beam  patterns  and  therefore  are  not  shown  for  signal  frequencies 
of  3  and  5  kHz.  A  sum  frequency  beam  pattern  for  a  signal  frequency  of 
9  kHz  is  shown  in  Fig.  6.9(c).  The  theoretical  result  from  Equation  ( 4 . 17" ) 
was  used  in  em  HP  9^30  basic  computer  program;  this  program  is  given  in 
Appendix  C.  Some  inconsistency  was  observed  in  the  experimental  begun 
patterns,  the  side  lobes  and  the  symmetry  being  somewhat  irregulsir  in  some 
cases.  This  irregularity  is  believed  to  be  due  to  the  presence  of  the 
support  structure  and  the  rotator.  Even  with  this  experimental  difficulty, 
the  agreement  between  the  theory  and  the  experiment  is  seen  to  be  good. 

Since  the  result  given  in  Equation  (4.l4)  is  considerably  simpler 
them  that  of  Equation  (4.17),  a  corapeirison  between  these  two  theoretical 
results  is  shown  in  Fig.  6.10  at  three  signal  frequencies  at  0=0.  The 
comparison  between  the  two  results  indicates  that  for  0=0  sufficient 
accuracy  is  achieved  by  using  the  simpler  equation  (Equation  4,l4).  The 
accuracy  of  Equation  (4.l4)  will  be  examined  further  when  the  parametric 
girray  beam  patterns  with  the  pump  transducer  misaligned  are  considered. 

The  second-order  sound  pressure  as  a  function  of  0'  for  the  sum 
and  difference  frequencies  with  a  >  kHz  signal  frequency  is  shown  in 
Fig.  6.11.  The  experimental  beam  pattern  was  plotted  carefully  on 
rectangular  chart  paper  with  an  angular  scale  of  3° /in.  to  ensure  adequate 
angular  resolution.  The  agreement  between  the  theory  and  experiment  was 
excellent.  Equation  (4.l4)  was  used  to  generate  the  theoretical  result. 

The  experimental  difference  frequency  beam  pattern  is  compared  to  the 
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FIGURE  6.9 

PARAMETRIC  RECEIVING  ARRAY  BEAM  PATTERNS 
FOR  A  NARROW  BEAM  PUMP  TRANSDUCER 


FIGURE  6.10 

COMPARISON  BETWEEN  THE  TWO  THEORETICAL  RESULTS 
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FIGURE  6.11 

BEAM  PATTERNS  FOR  THE  SUM  AND  DIFFERENCE  FREQUENCY 
SOUND  FIELDS  AS  PUMP  ANGLE  (S')  IS  VARIED 
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experimental  pump  beam  pattern  in  Fig.  6.12.  The  difference  frequency 
beam  pattern  has  the  same  beamwidth,  within  experimental  error,  as  the 
pump  beam  pattern. 

A  series  of  measurements  were  taken  with  the  pump  misaligned. 

The  pump  bearing  (0')  was  adjusted  to  the  3,  6,  or  10  dB  downpoint  on  the 
pump  beam  pattern.  The  sum  and  difference  frequency  parametric  array  beam 
patterns  for  a  signal  frequency  of  5  kHz  with  the  pump  misaligned  to  the 
-3  dB  point  (0'=1.2)  are  shown  in  Fig.  6.13.  The  sum  and  difference 
frequency  beam  patterns  are  no  longer  identical  or  symmetrical  as  they 
were  for  0  '=0,  nor  is  the  maximum  for  the  major  lobe  at  0=0.  One  side 
lobe  on  each  beam  pattern  was  higher  uhan  the  other  side  lobe.  The 
theoretical  beam  patterns  were  computed  using  Equation  (4.17).  The  sum 
and  difference  frequency  beam  patterns  for  the  pump  beam  at  -6  dB 
(0'=-1.6),  and  -10  dB  (0'=-2.O)  are  shown  in  Figs.  6.l4  and  6.15.  Again, 
the  two  theoretical  results  (Equations  4.l4  and  4.17)  are  compared  for 
0'=-1.6  in  Fig.  6.16.  In  each  of  these  beam  patterns  we  see  that  the 
angle  0,  for  which  the  array  response  is  a  maximum,  is  steered  toward  the 
direction  0'  is  rotated  for  the  difference  frequency  and  away  from  the 
direction  0'  is  rotated  for  the  sum  frequency.  This  property  of  the  sum 
and  difference  frequency  sound  fields  is  related  directly  to  the  idea  of 
"doppler"  angles.  These  results  indicate  that  the  sum  frequency  sound 
field  propagates  in  a  slightly  different  direction  than  the  original  pump 
sound  field,  whereas  the  difference  frequency  sound  field  propagates  in 
the  same  direction  as  the  original  pump  wave.  This  result  can  be 
demonstrated  more  dramatically  with  an  appropriate  coordinate  transforma¬ 
tion.  This  will  be  done  after  discussing  0  beam  patterns  with  6¥0. 


FIGURE  6.12 

COMPARISON  OF  THE  DIFFERENCE  FREQUENCY  SIDEBAND  BEAM  PATTERN 
WITH  THE  PUMP  BEAM  PATTERN  AS  THE  ANGLE  (S')  IS  VARIED 


(o)  SUM  FREQUENCY 


(b)  DIFFERENCE  FREQUENCY 


FIGURE  6.13 

PARAMETRIC  RECEIVING  ARRAY  BEAM  PATTERNS 
WITH  PUMP  MISALIGNED  ($' •  1.2°) 
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(a)  SUM  FREQUENCY 


(b)  DIFFERENCE  FREQUENCY 


FIGURE  6.14 

PARAMETRIC  RECEIVING  ARRAY  BEAM  PATTERNS 
WITH  PUMP  MISALIGNED  (0'=  -1.6°) 
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(a)  SUM  FREQUENCY  (b)  DIFFERENCE  FREQUENCY 


FIGURE  6.16 

COMPARISON  OF  THE  TWO  THEORETICAL  RESULTS  FOR 
PARAMETRIC  RECEIVING  ARRAY  BEAM  PATTERNS  FOR  0'=-1.6# 
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Beam  patterns  with  varying  0*  for  0=7°,  10°,  and  13°  for  a  5  kHz 
signal  frequency  are  shown  in  Figs.  6.17,  6. l8,  and  6.19.  The  value  for 
the  maximum  sum  and  difference  frequency  sound  pressures  is  seen  to  he 
slightly  different  than  zero.  For  the  sum  frequency,  the  maxima  were  at 
0 '=0.25°,  0.3°,  and  0.45°  for  0=7°,  10°,  and  13°,  respectively.  For  the 
difference  frequency  the  maxima  were  at  0'=-O.25,  -0.3,  and  -0.45  for 
0=7°,  10°,  and  13°,  respectively.  The  angle  0=13°  corresponded  to  the 
angle  for  the  first  side  lobe.  The  value  for  the  maximum  response  for 
0=0  is  marked  on  each  plot. 

Doppler  angles  using  the  geometry  shown  in  Fig.  6.20  will  now  be 
examined.  The  second-order  sound  field  will  be  observed  at  a  point  (L,0')« 
The  angle  of  the  low  frequency  plane  vector  with  respect  to  the  axiB  of  the 
pump  beam  pattern  is  0.  Plots  of  the  theoretical  sound  field  for  the  sum 
and  difference  frequency  components  as  a  function  of  0',  with  0=5°  and 
0=7°,  are  shown  in  Fig.  6.21.  The  signal  frequency  is  5  kHz  and  the  pump 
frequency  is  90  kHz.  The  sum  frequency  sound  wave  is  seen  to  propagate 
with  a  maximum  value  at  ('= 0,  as  for  0=5°  and  0'=O.4°  for  0=7°.  The 
difference  frequency  sound  wave  propagates  in  the  same  direction  as  the 
pump  wave. 

In  light  of  excellent  agreement  between  the  theoretical  and 

experimental  results,  there  is  little  doubt  that  any  significant  difference 

would  be  observed  for  an  experiment  done  using  the  coordinate  system  shown 

in  Fig.  6.20.  We  see  that  our  theory  and  experiment  is  consistent  with 

the  idea  of  interaction  occurring  at  nonzero  angles  between  the  direction 

of  propagation  of  the  two  waves.  Furthermore,  we  see  that  the  "doppler" 

15  49 

angles  predicted  by  Pridmore-Brown  and  Ingard,  Dean,  and  Lauvstad  and 
Tjotta,^0  can  indeed  be  observed. 
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FIGURE  6.17 

BEAM  PATTERNS  AT  THE  SIDEBAND  FREQUENCIES  FOR  0  =  7°  WITH  VARYING  0’ 
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FIGURE  6.19 

BEAM  PATTERNS  AT  THE  SIDEBAND  FREQUENCIES  FOR  9  -  13°  WITH  VARYING  O' 
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FIGURE  6.20 

\MINATI0N  OF  THE  DOPPLER  ANGLES 
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C .  Experiments  on  the  Parametric  Receiving  Array  vith  an  Omnidirectional 

Pump  and  a  Line  Receiver 

We  shall  now  reverse  the  high  frequency  transducers  and  examine 
the  properties  of  a  parametric  receiving  array  with  a  narrowbeam  receiver. 
We  first  examine  sum  frequency  beam  patterns  for  the  parametric  receiving 
array  at  0  '=0.  The  beam  patterns  with  signal  frequencies  of  5  and  9  kHz 
are  shown  in  Fig.  6.22.  !Diese  beam  patterns  are  very  Bimilar  to  those 
for  a  narrowbeara  pump.  Equation  (4.32)  was  used  for  calculating  the 
theoretical  beam  patterns.  The  agreement  with  the  theoretical  results  was 
excellent.  The  support  for  the  omnidirectional  transducer  was  a  thin  rod 
which  did  little  to  affect  the  low  frequency  sound  field. 

Difference  frequency  beam  patterns  with  the  receiver  misaligned 
to  the  3  <JB  (0'=±1.3°)  and  6  dB  (0*=±1.7°)  points  at  a  signal  frequency  of 
5  kHz  are  shown  in  Figs.  6.23  and  6.24,  respectively.  These  beam  patterns 
resemble  those  found  for  the  narrowbeam  pump.  The  theoretical  beam 
patterns  were  generated  using  Equation  (4.32).  The  array  beam  pattern  at 
the  sum  frequency  for  the  narrowbeam  receiver  at  -3  dB  (0'=l.l)  and  at 
-6  dB  (0'=1.53)  are  shown  in  Fig.  6.25* 

Beam  patterns  obtained  by  varying  0'  with  0=0  are  shown  for  the 
6um  and  difference  frequencies  in  Fig.  6.26.  There  is  an  important 
difference  in  these  beam  patterns  from  the  original  pump  beam  patterns. 

'The  beamwidth  of  these  patterns  were  found  to  be  that  of  the  receiver  at 
the  sum  and  difference  frequencies  rather  than  at  the  pump  frequency,  as 
was  observed  with  the  narrowbeam  pump.  The  theoretical  results  were 
plotted  using  Equation  (4.29). 

Receiver  beam  patterns  with  varying  0'  with  an  array  angle  of 
0=7°  and  13°  are  shown  in  Figs.  6.27  and  6.28,  respectively.  Again,  a 
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(o)  5  kHz  SUM  FREQUENCY 


(b)  9  kHz  SUM  FREQUENCY 


FIGURE  6.22 

SUM  FREQUENCY  PARAMETRIC  RECEIVING  ARRAY  BEAM  PATTERNS 


1.3°  DIFFERENCE  FREQUENCY  (b)  B'  -  -1.3°  DIFFERENCE  FREQUENCY 


FIGURE  6.23 

DIFFERENCE  FREQUENCY  ARRAY  BEAM  PATTERNS 
WITH  THE  RECEIVER  MISALIGNED  BY  9' 
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FIGURE  6.24 

DIFFERENCE  FREQUENCY  ARRAY  BEAM  PATTERNS 
WITH  THE  RECEIVER  MISALIGNED  BY  9' 
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FIGURE  6.25 

SUM  FREQUENCY  ARRAY  BEAM  PATTERNS 
WITH  THE  RECEIVER  MISALIGNED  BY  9' 
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FIGURE  6.27 

SIDEBAND  RECEIVER  BEAM  PATTERNS  FOR  0  *  V 
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FIGURE  6.28 

SIDEBAND  RECEIVER  BEAM  PATTERNS  FOR  8  ■  13* 
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slight  change  in  the  angle  9 '  for  the  maximum  was  slightly  different 
than  zero. 

D.  Experiments  on  the  Parametric  Receiving  Array  with  a  Small  Pump 
Transducer 

A  small  transducer  with  a  measured  beam  pattern  at  $0  kHz  as 
shown  in  Fig.  6.29  was  used  for  the  pump  transducer.  The  beam  pattern  for 
the  transducer  was  found  to  be  very  similar  to  that  of  a  4  in.  x  4  in. 
square  aperture  used  for  the  theoretical  result  in  Fig.  6.29.  A  beam 
pattern  for  the  parametric  array  with  a  signal  frequency  of  5  kHz  is 
shown  in  Fig.  6. 30.  As  with  the  narrowbeam  pump  transducer,  some  irregu¬ 
larity  in  the  beam  pattern  was  observed.  This  irregularity  is  attributed 
to  the  presence  of  the  rotator  and  mounting  structure.  Equation  (4.29) 
was  used  to  generate  the  theoretical  beam  pattern.  The  sum  frequency  beam 
pattern  obtained  by  varying  6 '  with  the  array  aligned  with  the  signal  wave 
is  shown  in  Fig.  6. 31.  Again,  Equation  (4.29)  was  used  to  generate  the 
theoretical  result.  The  beam  pattern  closely  resembles  that  of  the  pump 
transducer  operated  at  90  kHz  as  is  predicted  by  the  theory. 
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FIGURE  6.30 

SUM  FREQUENCY  PARAMETRIC  ARRAY  BEAM  PATTERN 
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FIGURE  6.31 

SUM  FREQUENCY  BEAM  PATTERN  WITH  VARYING  6' 


VII.  SUMMARY  AND  CONCLUSIONS 


Two  objectives  were  undertaken  in  this  study  of  the  parametric 
receiving  array.  They  were  (l)  to  relate  the  parametric  receiving  array 
to  the  problem  of  interaction  of  sound  by  sound  at  nonzero  angles,  and 
(2)  to  find  a  solution  for  the  second-order  sound  field  for  a  parametric 
receiving  array  with  pump  sources  which  ranged  from  omnidirectional  to 
very  narrowbeam  acoustic  generators.  Each  of  these  objectives  was 
accomplished  using  the  perturbation  solution  developed  by  Westervelt.  The 
second-order  sound  pressure  was  obtained  by  solving  an  inhomogeneous 
wave  equation  with  a  source  function  expressed  in  terms  of  a  product  of 
first-order  sound  field  variables.  In  the  derivation  of  the  wave  equation, 
the  effects  of  viscosity  and  thermal  conductivity  were  ignored.  The 
effects  of  absorption  were  included  in  the  expression  of  the  first-order 
field  variables. 

Thu  interaction  of  two  plane  waves  over  all  space  in  a  lossless 
medium  was  considered  in  Chapter  III.  Westervelt 's  1957  result  was 
reproduced.  This  result  predicted  an  infinite  sound  pressure  when  the 
two  plane  waves  are  propagating  in  the  same  direction.  When  the  sound 
waves  are  not  propagating  in  the  same  directions,  the  second-order  sound 
field  was  found  to  be  a  function  of  the  local  value  for  the  first-order 
sound  field.  A  more  reasonable  solution  was  then  obtained  by  placing  a 
boundary  condition  on  one  of  the  sound  fields.  For  the  parametric 
receiving  array,  the  sound  wave  with  the  boundary  condition  was  assumed 
to  have  a  frequency  f^,  such  that  f^  was  considerably  greater  than 
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f2,  the  signal  frequency.  This  case  represented  a  parametric  receiving 
array  with  an  infinite  planar  piston  source  for  a  pump.  The  second-order 
sound  pressure  was  found  to  be  proportional  to  the  distance  of  the  observer 
from  the  planar  boundary  and  hr'l  the  directivity  function  of  a  truncated 
end-fire  array.  Berktay  observed  a  similar  result  for  the  parametric 
receiving  array  with  the  receiving  transducer  in  the  nearfield  of  the  pump 
transducer.  Ihe  same  end-fire  function  was  predicted  for  an  omnidirec¬ 
tional  pump  transducer  in  our  study.  The  interaction  of  sound  with  sound 
was  also  considered  in  Chapter  VI  when  we  observed  the  properties  of  the 
second-order  sound  field  at  nonzero  angles  for  the  array.  A  modified  form 
for  the  doppler  angles  was  predicted  by  the  theory  and  observed  by  the 
experiment.  The  form  for  the  doppler  angles  was  modified  because  the  high 
frequency  pump  wave  propagated  as  a  spherical  rather  than  a  planar  wave. 

In  the  present  study,  the  observer  waB  always  in  the  region  of  interaction. 
It  can  be  concluded  that  cumulative  Interaction  of  sound  by  sound  at 
nonzero  angles  does  occur  within  the  zone  of  interaction.  No  conclusion 
about  the  presence  of  a  second-order  sound  wave  outside  the  region  of 
interaction  can  be  drawn  since,  strictly  speaking,  the  observer  waB  always 
in  the  zone  of  interaction. 

An  interesting  special  problem,  the  parametric  receiving  array 
with  a  point  source  pump,  was  solved  in  Chapter  III.  This  solution  is  of 
particular  significance  for  two  reasons  (l)  this  case  represents  a 
situation  which  can  be  modeled  both  theoretically  and  experimentally,  and 
(2)  the  solution  can  be  used  to  solve  parametric  receiving  array  problems 
with  more  complex  pump  source  distributions.  In  our  analysis,  it  has  been 
assumed  that  the  pump  wave  and  signal  wave  amplitudes  are  small  enough  to 
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use  a  small-signal  nonlinear  acoustics  or  quasilinear  solution  to  the 
problem.  The  point  source  solution  was  used  as  a  starting  point  to 
obtain  the  solution  for  a  line  source  pump.  1Mb  summation  vas  legitimate 
only  because  the  amplitude  of  the  pump  wave  vas  assumed  to  be  small  enough 
to  preclude  finite  amplitude  attenuation  of  the  pump  vave.  The  effects  of 
absorption  were  included  in  the  description  of  the  firec-order  sound  field. 
The  solution  for  the  point  source  pump  wave  was  obtained  using  a  two- 
dimensional  stationary  phase  integral  evaluation  technique.  This  approach 
represents  an  interesting  new  means  to  solve  three-dimensional  nonlinear 
acoustics  problems.  To  the  author's  knowledge,  this  application  of  the 
two-dimensional  integral  technique  is  the  first  for  a  nonlinear  acoustics 
problem.  The  parametric  receiving  array  with  a  point  source  pump  was 
found  to  have  a  directivity  function  equal  to  that  of  an  equivalent 
end-fire  array.  When  the  pump  source  was  misaligned  with  respect  to  the 
receiver,  the  directivity  function  was  dramatically  changed.  The  beam 
patterns  found  by  rotating  the  pump  with  the  array  aligned  with  the  low 
frequency  sound  field  was  found  to  be  almost  identical  to  that  of  the  pump 
operating  at  the  pump  frequency.  These  properties  of  the  parametric 
receiving  array  with  a  line  source  pump  were  observed  both  theoretically 
and  experimentally. 

The  line  transducer  was  used  as  a  receiver  with  an  omnidirectional 
pump  source.  Again,  the  array  beam  patterns  were  found  to  be  narrower 
with  lower  side  lobes  than  the  equivalent  end-fire  array.  Beam  patterns 
observed  by  rotating  the  receiver  transducer  were  found  to  be  the  same  as 
beam  patterns  made  at  the  sum  or  difference  frequency.  Again,  the 
theoretical  and  experimental  results  agreed  closely. 
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When  a  transducer  with  dimensions  small  compared  to  a  signal 
wavelength  was  used  as  a  pump  transducer,  the  properties  of  the  parametric 

receiving  array  were  only  slightly  different  from  those  predicted  by 

32  33  39  4o 

Berktay  and  A1  Temimi,  *  Berktay  and  Shooter,  7  and  Rogers  et  al. 

In  this  study  of  the  parametric  receiving  array,  only  the 
cumulative  contribution  to  the  second-order  pressure  was  considered.  The 
second-order  effects  at  the  pump  and  at  the  receiver  were  ignored.  Since 
the  parametric  receiver  offers  potential  usefulness  as  a  low  frequency 
narrowbeam  receiver,  additional  analysis  of  these  local  pressure  terms  is 
recommended  since  these  terms  could  reduce  the  directivity  of  the  parametric 
receiving  array.  The  effects  of  noise  at  the  side  band  frequency  has  not 
been  considered  in  the  present  study.  Side  band  noise  will  be  generated 
by  the  electronics  associated  with  both  the  pump  and  receiver.  Further¬ 
more,  side  band  noise  will  be  generated  if  the  direct  or  reflected  acoustic 
path  length  from  the  pump  to  the  receiver  varies  with  time.  The  amplitude 
of  this  side  band  noise  needs  to  be  determined  under  the  actual  conditions 
in  which  the  parametric  receiving  array  is  considered  for  use. 

The  effects  of  irhomogeneities  on  the  water  medium  likewise 
needs  to  be  considered.  Inhomogeneities  in  the  region  between  the  pump 
and  receiver  could  reduce  the  amplitude  of  the  desired  side  band  pressure 
and  also  degrade  the  beam  pattern  of  the  parametric  receiving  array. 

It  is  hoped  that  this  study  has  laid  a  firm  foundation  for  the 
understanding  of  the  parametric  receiving  array  with  various  pump  and 
receiver  transducers.  Likewise,  it  is  hoped  that  this  work  has  shed  some 
light  on  the  problem  of  sound  scattered  by  sound. 


APPENDIX  A 

TWO-DIMENSIONAL  STATIONARY  PHASE  SOLUTION  OF  THE  SCATTERING  INTEGRAL 


In  Section  D  of  Chapter  III,  we  had  the  Green's  function  integral 


solution 


P  = 
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Jtflp  c 
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P11P12 


+  z 


(A.l) 


•  expj-a^x2  +  y2  +  z2  -  a2(xo  cos  Q  +  y  sin  0)  -  a±>/ (xq-x)2  +  y2  +  £ 

•  exp^C^^x2  +  y2  +  z2  i  1^  x  cos  6  ±  y  sin  0)j 

.  expJjCk^kg)  >/(xQ-x)2  +  +  z2  -  ^(a^ta^Jtj  dxdydz 

For  a  given  value  of  x,  the  integrand  amplitude  term  is  a  slowly  varying 
function  of  y  and  z.  On  the  other  hand,  the  phase  varies  rapidly  with  y 
and  z.  These  properties  of  the  integrand  suggest  the  use  of  the  method  of 
stationary  phase. 

Double  integrals  of  the  fo^n^,^,^ 

ff  S(y,2)eJllf(y’2)  iyiz  (A.2) 

can  be  evaluated  by  the  method  of  stationary  phase  provided  S(y,z)  is  a 
slowly  varying  function  of  y  and  z,  and  kf(y,z)  is  a  rapidly  varying 
function  of  y  and  z.  The  contributions  to  the  asymtotic  expansion  of  the 
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integral  come  only  from  regions  in  the  vicinity  of  certain  critical  points. 
There  are  three  types  of  critical  points.  The  critical  points  of  interest 
for  our  problem  are  those  within  the  domain  of  integration  at  which 


df  _  df  0 


(A.  3) 


Then  near  the  critical  point  (yo,zQ)  we  can  expand  f(y,z)  in  a  Taylor 
series  such  that  we  have 

f(y,z)  a  f(y0,zQ)  +  \  a(y-y0)2  +  ^  £(z-zo)2  +  7(y-yQ)  (z-zq)  ,  (A. 4) 


where 


a  = 


it 

*y2 


P~dz2 


<£f 

7  =  3y¥ 


the  partial  derivatives  all  being  evaluated  at  (yo,zQ).  We  can  now  choose 


new  variables  of  integration  such  that 


f(y»z)  =  f(y0>z0)  +  £a  2  +  z  V2  +  rSn 


(A.5) 


The  asynrtotic  approximation  to  Equation  (A.l)  is 


S(yo,zo)eikf(yo^zo)  j  expj^  ilcCa2  +^T]2-^rnl )J  d|dTj  = 


ikf(y  ,z  ) 

2nj  a  S(y  ,z  )e 
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where  the  positive  root  is  taken  and  a=+l,  -1,  or  -j  according  to  whether 

ap>y2fd>0  op>)P,a <0  ;  (A. 7) 


or  a {J<7  . 
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Then  -he  functions  f(y,z)  and  S(y,z)  are  given  by 

kl  ./T~  2  .  2  .  V 

rr  vx  +  y  +  z  ±  r — r-; 


f(y>z)  ^ 


ki  1  *2 


n  rs-;  2  .  2 

v(x  -x) 


^±1^  8lne+  v  1  vxr  +  y  +z" 


and 


j(<0  ±0i,)2  (1+B/2A)  PuPjo 
S(y,z)  =  - - - y 


Wlp  c 
0  o 


(A. 8) 


exp|-a2(x  cos  6  +  y  sin  0)  -  a^J x^  +  +  z^ 


-  cr 


>/(xq-x)2  +  y2  +  z2  ±  kgX  cos  0 


It  can  be  seen  that  f(y,z)  is  a  rapidly  varying  function  of  y  and  z  while 
S(y,z)  is  a  slowly  varying  function  of  y  and  z. 

To  find  the  critical  point,  we  must  find  the  point  at  which 

df  _  8f 
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For  the  z  derivative  we  have 
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(A. 9) 


Jt  "  \g  '"'S'  ^  5 

v(x  -x)  +  y  +  z 


=  0  if  z  =  0 


and  for  the  y  derivative  we  have 


kl±k2 


sin  6  + 


7^7777 


=  0  .  (A.  10) 
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Solving  for  y  we  have 


y  * 


*2 


sin  9 


*1  ,  1 
(kjikgjx  (xQ-x) 


kg  sin  6  (x)  (xQ-x) 
T  k1(xQ-x)  +  xCk^kg) 


kg  sin  9  (x)  (xQ-x) 

(A.  11) 


We  now  can  obtain  the  values  of  the  second  derivatives  of  f  evaluated  at 
the  critical  point.  We  have 
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(A. 12) 

k^2 

y2 

(kjtkg)  (x2+y2+z2)5/2 

((V*)2+y2+z2)5/2 

(A.  13) 
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Also  we  have 
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The  function  f  evaluated  at  the  critical  point  is 


f<y0>z0) 


kl  ±  *2 


'x  +  yo  ±  y  8in  6  +  v(x-xo)  +  y 


.  .2  2  ,  2  .2  ,2.  (A.15) 
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The  second-order  pressure  is  then 
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APPENDIX  B 
TRANSDUCERS 

The  five  transducers  used  for  the  parametric  receiving  array 
experiments  are  described  in  the  following  paragraphs. 

A.  Rectangular  Transducer 

This  transducer  was  constructed  using  elements  for  use  as  a 
narrowbeam  receiver.  In  our  experiments,  the  transducer  was  used  as  both 
a  projector  and  as  a  hydrophone.  The  dimensions  of  the  active  face  were 
17.5  in.  by  2.75  in.  The  end  elements  are  reduced  in  amplitude  to  provide 
amplitude  shading  to  reduce  the  side  lobes.  Beam  patterns  for  this 
transducer  at  90  kHz  are  shown  in  Fig.  6.4. 

B.  Small  Standard  Transducer 

This  is  a  U.  S.  Navy  standard  TR-129  transducer  designed  as  a 
small  general  purpose  transducer.  It  was  constructed  of  small  ceramic 
cylinders  designed  to  radiate  omnidirectionally  in  the  plane  perpendicular 
to  the  axis  of  the  cylinders.  The  cylinders  have  a  circumferential 
resonance  at  90  kHz. 

C.  Standard  Hydrophone 

This  is  a  U.  S.  Navy  standard  TR-225  transducer  developed  by  the 
Naval  Research  Laboratory,  Underwater  Sound  Reference  Division,  Orlando, 
Florida.  It  is  constructed  with  lead-zirconate/ lead-titanate  cylinders. 
Like  the  previous  transducer,  it  is  designed  to  be  omnidirectional  in  the 
plane  perpendicular  to  the  axis  of  the  cylinders.  The  transducer  has  a 
nominal  freefield  open  circuit  voltage  sensitivity  of  -203  dB  re  1  v/l 

and  has  a  flat  response  (±1  dB)  over  the  frequency  range  from  1  to  20  kHz. 
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D.  Transducer 

This  transducer  was  used  as  a  pump  and  as  a  receiver.  Ihls 
transducer  was  constructed  at  the  Applied  Research  Laboratories  and  has 
been  assigned  the  number  DRL-123.  It  was  constructed  as  a  mosaic  of 
Channelite  5*+00  ceramic  cubes.  The  transducer  had  a  resonant  frequency 
of  90  kHz. 

E.  Hydrophone 

This  is  a  U.  S.  Navy  standard  H-23  hydrophone  developed  and 
supplied  by  the  Naval  Research  Laboratory,  Underwater  Sound  Reference 
Division,  Orlando,  Florida.  It  was  constructed  using  eight  lithium 
sulfate  crystals.  It  has  a  follower  amplifier  to  drive  the  cable.  It  haB 
a  nominal  sensitivity  of  -l88  dB  re  1  V/l  uPa  at  90  kHz. 


APPENDIX  C 
COMPUTER  PROGRAM 


A  computer  program  was  written  by  J.  Kodosky  and  W.  C.  Nowlin 
for  use  on  the  Hewlett-Packard  9^30  Programmable  Calculator  in  BASIC 
language.  The  program  was  designed  to  plot  0  or  0'  beam  patterns  using 
Equations  (4.l4),  (4.17),  ('•-.29),  (4.32).  The  choice  of  equation 

calculated  was  input  selectable.  The  program  allowed  rotation  of  the  zero 
for  0  and  0’  for  alignment  of  the  theoretical  beam  patterns  with  the 
experimental  beam  pattern.  The  amplitude  of  the  theoretical  beam  pattern 
could  also  be  adjusted  to  account  for  amplitude  variations  in  the  experi¬ 


mental  results. 
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